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Programmgesteuette. digitale Rechengerate 
_ (elektronische Rechenmaschinen) 


Zusammenfassender Bericht (II. Teil) 
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§ 3. Arithmetische Prinzipien 


3.1. Zahlsysteme 


Unter einem Zahlsystem versteht man allgemein die Art und Weise, eine — 
positiv reelle Zahl durch Zeichen (Ziffern) darzustellen, d.h. zu verziffern.® 
Dabei ist das Dezimalsystem nur ein Sonderfall (6 =10) der allgemeineren 
Darstellung 


M 
x= 3" x, Be, (3.1) 


wobei die ganze Zahl B = 2 die Basis des Zahlsystems heiBt und die ebenfalls 
ganzen Zahlen x, (0 < x, < B) die Ziffern der Zahl x sind. 

Beim praktischen Rechnen kénnen immer nur endliche B-albriiche verwen- 
det werden, die dann von den darzustellenden Zahlen « um gewisse Betrage 
6x abweichen; deshalb interessiert die mit einer gewissen Stellenzahl N erreichte 
relative Genauigkeit y = x/6x. Da durch geeignete Abrundung immer erreicht 
werden kann, daB 6x héchstens eine halbe Einheit der letzten Stelle ist und 
eine N-stellige B-alzahl stets mindestens B’~1 Einheiten der letzten Stelle 
betragt, gilt 2 2 BY}. 

Um also auch in ungiinstigen Fallen eine relative Genauigkeit yp) zu sichern, 
sind 


log y,/2 
errs en (3.2) 
wesentliche Stellen notwendig. 

Wahrenddem nun bisher die Zahl zehn als Basis fiir die Verzifferung unbe- 
stritten war, kommen fiir den Bau programmgesteuerter Rechenmaschinen 
doch auch andere Zahlsysteme in Betracht. Zur Verzifferung mit der Basis B 
braucht es B verschiedene Ziffern (0,1,..., B — 1), die in einer Rechenma- 
schine durch B verschiedene Zustande eines geeigneten mechanischen oder 
elektrischen Systems dargestellt werden kénnen. So haben die Zahlrader in den 
Burorechenmaschinen 10 ausgezeichnete Stellungen, die den Ziffern 0 bis 9 
zugeordnet sind. 

Fiir eine schnelle Maschine kommen hingegen nur elektrische Schaltelemente, 
zum Beispiel Relais und Elektronenréhren, in Betracht, und diese haben nur 
zwei stabile Zusténde. Somit ist fiir solche Maschinen das Zahlsystem mit der 
Basis 6b = 2, das sogenannte Dualsystem, bevorzugt, welches tatsachlich bei 
einigen fertiggestellten oder im Bau befindlichen Maschinen zur Anwendung 
kommt?) : 

Da das Dualsystem eine rein interne Angelegenheit dieser Maschinen ist — 
man denkt wohl kaum daran, auch im taglichen Leben das Dezimalsystem zu 


1) Vgl. die Zusammenstellung der Maschinen am Schlu8 von § 2, S. 291-294, 


# a 


x. 


a 


Ss 
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verdrangen — miussen die Ausgangswerte und Konstanten eines Problems 
zuerst ins Dualsystem tibertragen werden, und am Ende der Rechnung hat die 
Riickverwandlung der Resultate ins Dezimalsystem zu erfolgen (vgl. § 3.7). 


Die Vorziige des Dualsystems sind durch die Tatsache begriindet, da8 man 
mit den beiden Ziffern 0 und L1) auskommt, jede Ziffer einer Zahl ist damit 
immer entweder 0 oder L. Diese Alternative entspricht nicht nur den elektrischen 
Schaltungen vortrefflich (geéffnetes oder geschlossenes Relais), sondern stellt 

auch die Briicke zum Aussagenkalkiil her?) (richtig oder falsch). 

Eine fiir die Praxis wichtige Eigenschaft ergibt sich aus Formel (3.1), vermége 
der sich eine Zahl aus ihren Ziffern aufbaut: Nur beim Dualsystem ist ¥ immer 
eine Summe der reinen Potenzen der Basis, z. B. LOLOL = 24+ 22+ 29= 21; bei 

allen andern Zahlsystemen sind die Faktoren x, im allgemeinen gréBer als 1. 

Ferner zeigt es sich, daB das Dualsystem im Mittel die kleinsten Quersummen 
ergibt, was fiir Multiplikation und Division von Bedeutung ist, denn die Dauer 
dieser Operationen ist durch die Quersumme des Multiplikators bzw. Quotienten 
bestimmt. Weil jede Stelle eines B-albruches im Mittel (B— 1)/2 zur Quersumme 

-beitragt und die zur Erreichung einer relativen Genauigkeit y, notwendige Stellen- 
zahl nach (3.2) im wesentlichen log y,/log B ist, folgt, daB die mittlere zu er- 
wartende Quersumme bei festem, vorgeschriebenem yj zu (B— 1)/log B propor- 
tional ist. Dieser Ausdruck ist fiir B = 2 am kleinsten, namlich zirka 1,44, wahrend 

‘sich fiir das Dezimalsystem zirka 3,9 ergibt. 

Zu den Nachteilen des Dualsystems mu8 einmal die Notwendigkeit des Uber- 
setzens gerechnet werden, ferner ist das Ziffernbild einer Zahl wenig differenziert, 
weil nur die Ziffern 0 und L auftreten. Es ist deshalb nicht leicht, eine Zahl mit 
einem Blick zu erfassen, um so mehr, als im Dualsystem wesentlich mehr Stellen 
notwendig sind, um die gleiche Genauigkeit zu gewdhrleisten (einer 12stelligen 
Dezimalzahl entspricht eine 40stellige Dualzahl). 


3.2. Duale Verschliisselung 


Da die Verwendung des reinen Dualsystems in einer Rechenmaschine viele 
Nachteile mit sich bringt, haben namhafte Konstrukteure am Dezimalsystem 
festgehalten. Um aber dennoch von den Vorteilen des Dualsystems profitieren zu 
kénnen, stellen sie wenigstens die einzelnen Ziffern einer Zahl durch Dualzahlen 
dar. Es wird zu diesem Zweck eine Vorschrift eingefiihrt, welche jeder ganzen 
Zahl x (0 < x<10) eindeutig eine ebenfalls ganze Zahl ¢(x) [0 S#(x) < 16] 
zuordnet. Diese wird als vierstellige Dualzahl geschrieben und die der Ziffer x 
zugeordnete Tetrade genannt. 

Die nachstliegende aller méglichen Vorschriften ist durch die Funktion 
t(x) = x bestimmt (sogenannte direkte Verschliisselung)*). Bei dieser wird also 
einfach jede Ziffer ins Dualsystem tibersetzt, z. B. 


6295 > OLLO’00L0’LOOL’0LOL. 


1) Die Dual-Eins wird zur Unterscheidung von der Dezima!-Eins mit L bezeichnet. 


2) Vgl. C. E. SHANNON [48]. 
3) Angewendet bei Mark II [3] und SSEC (IBM) [82]. 


ee 
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Ferner wurde beim «Complex Computer»!), welcher allerdings keine eigent-_ 
liche programmgesteuerte Rechenmaschine ist, die sogenannte Dreiexzep-Ver-_ 
schliisselung verwendet: ¢(x) = x + 3, oder 


BT 


0 00LL, Le 0L00.. 25- OLOL i. Aas a Eee 


Diese zeichnet sich unter anderem dadurch aus, daB8 niemals alle Dualzif- 
fern einer verschliisselten Zahl 0 oder alle L sein kénnen, so daB sich das Ver- 
sagen von elektrischen Schaltelementen sofort anzeigt. 

Wir wollen uns nun sofort dem allgemeinen Fall zuwenden: Die Gesamtzahl 
aller méglichen Zuordnungen «-—> d(x) betragt offenbar 161°, worin allerdings 
auch die nicht eindeutig umkehrbaren inbegriffen sind. Um die Ubersicht nicht 
zu verlieren, scheidet man aus dieser Menge durch geeignete Postulate die- 
jenigen Funktionen ¢(«) aus, die zu unzweckmaBigen Verschliisselungen fihren 
wiirden. Diese Postulate sind nach AIKEN®): 


A. Verschiedenen Ziffern sollen verschiedene Tetraden zugeordnet werden, d. h. 
t(x) soll jeden Wert héchstens einmal annehmen. 


B. Die gréBere Ziffer soll auch durch die groBere Dualzahl dargestellt werden, 
d.h. aus x < y soll ¢(x) < ¢(y) folgen. 


C. Wenn sich zwei Ziffern x und y auf 9 erganzen, sollen auch die zugeordneten 
Tetraden komplementar sein, also durch Vertauschung von 0 und L aus- 
einander hervorgehen, d. h. es soll dann ¢(x) + ¢(y) = 15 gelten. 


Aus diesen Postulaten, denen noch 56 Verschliisselungen geniigen, folgt 
librigens, da8B sich die Ziffern 0 bis 4 einerseits und 5 bis 9 anderseits durch die 
erste Dualziffer der zugeordneten Tetraden unterscheiden; diese ist namlich 
Ove) <8) tine << 5.und Lid 8) fires = 5. 

D. Die geraden (bzw. ungeraden) Ziffern sollen durch gerade (bzw. ungerade) 

Tetraden (oder umgekehrt) dargestellt werden, so daB die letzte Dualstelle 

einer Tetrade die Paritat der dargestellten Ziffer bestimmt. 


E. Es soll méglich sein, jeder der vier Dualstellen einer Tetrade geeignete 
Gewichte zuzuordnen, so da die gewogene Quersumme der Tetrade die dar- 
gestellte Ziffer ergibt; es soll sich also x auf diese einfache Weise aus f(x) 
berechnen lassen. (Zum Beispiel haben die vier Dualstellen bei der anfangs 
erwahnten direkten Verschliisselung die Gewichte 8, 4, 2, 1.) 


Die Zahl der verschiedenen Verschliisselungen, welche diesen fiinf Axiomen 
gentigen (die direkte Verschliisselung verletzt das Postulat C und die Drei- 
exzeBverschliisselung geniigt der Bedingung E nicht), betragt noch 4, unter 
denen sich die Aikensche’ oder 2-4-2-1-Verschliisselung durch besondere Ein- 
fachheit auszeichnet: 


1) Beschreibung dieser Maschine: Bell Laboratories Record TONGS 25S: 5/6 (Oktober 1940). 
2) Persénliche Mitteilung. 


. 
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og es fir wa 
elo firs 2 25 | 
oder: 
0+ 0000 5> LOLL 
1+ 000L 6+>LLOO (3:3) 
2> 00LO 7> LLOL 
3> OOLL 8—> LLLO 
4—> 0LO0 oS ELL | 


Die Gewichte sind hier 2, 4, 2, 1; zum Beispiel ist also7 (LLOL) = 2+4+0+41. 

Rechnen mit dual verschliisselten Zahlen. Dariiber sagen die fiinf Axiome 
gar nichts aus. Das Rechnen kann nun immer auf die Addition zuriickgefihrt 
werden!), welche ziffernweise durchgefiihrt wird, so daB man sich darauf be- 
schranken kann, die Addition von einzelnen Ziffern in Betracht zu ziehen. Da 
diese aber Ausschnitte aus vielstelligen Zahlen sind, miissen noch die eventuell 
erfolgenden Zehnertibertrage von rechts her und nach links beriicksichtigt 
werden. ; 

Die Aufgabe, die Summe zweier Ziffern x und y zu bilden, besteht also 
darin, aus den Tetraden ¢(x) und ¢(y) und dem eventuell von rechts kommen- 
den Ubertrag u (es ist « = 1 oder 0, je nachdem ein solcher Ubertrag erfolgt 
oder nicht) die folgende Tetrade zu bilden: 


tix+ y+ u) falls x + y+u< 10 
(3.4) 
16+¢*+v+u—10) falls x+y+u2Z10 


Der Summand 16 im zweiten Fall entspricht dem Zehneriibertrag, der zur 
nachst hoheren Tetrade zu erfolgen hat und der, dual aufgefaBt, einen Wert von 
16 Einheiten der betrachteten Tetrade besitzt (vgl. Fig. 2). 


Tetrade 10"t! Tetrade 10” 
Lg Seyehoobon teva Ge Sale xX X XX x X X X 
PeSUMINANC nice ios ONO x X X X 
ISOMmeK LUNE )s ee. x XXX xX XX 
Summe ohne Ubertrag OOo LL Eh LOL Liz. B:) 

(oe 

Summe mit Ubertrag . OLO0 Leon ie 

Fig. 2 

Zehneriibertrag 


1) Vel. § 3.5. 
2) Vel. weiter unten. 


Man geht nun meistens so vor, da8 man zunachst in einem dualen Addier- 
werk die Summe ¢(x) + ¢(y) + u bildet; falls dies nicht das gema8 (3.4) gefor- 
derte Resultat ergibt, hat noch eine Korrektur zu erfolgen, was je nach der ge- 
wahlten Verschliisselung mehr oder weniger kompliziert ist (vgl. auch § 5.11). — 
Um diese Korrekturen zu bestimmen, vergleichen wir die Summe ¢(x) + ¢(y) + u_ 
(Haben) mit den Sollbetragen von (3.4), wobei fiir ¢(x) die speziellen Funk-_— 
tionen eingesetzt werden: ; 


i 
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a) Direkte Verschliisselung 


Soll: x+y+u fir *«+y+tu< 10 


Haben: pty, 
64+x“+y+u4 fir *+y+u210 Os ue anaes 


Das heiBt: Die Summe ¢(x) + ¢(y) + wu ist richtig fiir «+ y+ u<10, andern- 
falls hat man die Korrektur 6 = OLLO zu addieren (worauf dann ein Ubertrag 
zur nachst hdheren Tetrade erfolgt). 


b) Dretexze6-Verschliisselung 


Soll: x+y+u+3 fir x*«+y+u<10 


xt+ty+u+9 fiir x+y +u 210 Haben : x+ytu+6. 


Die Korrektur ist also — 3 im ersten, + 3 im zweiten Fall. Diese beiden Falle 
unterscheiden sich offenbar dadurch, da8 der Habenbetrag im ersten Fall < 16, 
im zweiten = 16 ist, d. h. daB genau im zweiten Fall schon ohne die Korrektur 
ein Ubertrag zur nachst hdheren Tetrade erfolgt, welcher richtig ist und durch 
die Korrektur nicht mehr beeinflu8t wird. 


c) Atken-Verschliisselung 


Soll: Haben: 

mt au He Pa eS, Dade ar fiir x,y <5 

te Y-Pw 6 fir 5S See ye 15) 4 ye our = 5, Ves 
i He Sy ae 

Ary ut 12 fir I Set yt uw, e+ytuti2 fir x,y S5. 


Korrekturen sind also in folgenden Fallen notwendig: 


Be Sa a ap = 5: Korrektur + 6 , 
*%Y25, x«+y+u<15: Korrektur—6. 


Interessanterweise ist die Summe ¢(x) + ¢(y) + « gerade dann korrekturbe- 
dirftig, wenn sie, von einem eventuellen Ubertrag zur nachst hoheren Stelle 
abgesehen, eine Pseudodezimale ist, d. h. eine der 6 nicht zur Darstellung der 
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10 Ziffern verwendeten vierstelligen Dualzahlen: 


2 == ULOL 27 == 1000) 
6* = QOLLO 3* = LOOL 
W lag — ia 0 feet O10 


Wenn also die Summe #(x) + ¢(y) + wu eine richtige Tetrade ist, so bedarf sie 
keiner Korrektur mehr. 


Beispiel: 
= 5> LOLk IO l 
6—> LLOO 8+ LLLO 
L 0LLL-> Pseudodezimale 7* 15<L LOLL = richtige Dezimale 5 
+ Ubertrag, also 4+ Ubertrag, keine 
— 0LLO Korrektur — 0LLO Korrektur 


11<L0O00L 


Eine besondere Art der Darstellung der Ziffern (Biquinary Notation), welche 
den roémischen Ziffern gleicht, wurde im Bell Computer!) verwendet. 


3.3. Darstellung der Zahlen in der Maschine*) 


__ Im folgenden soll davon abgesehen werden, da die zehn Ziffern des Dezi- 
malsystems durch Tetraden dargestellt werden und da nach jeder Addition 
Korrekturen notwendig sind. 

In jeder digitalen Rechenmaschine kénnen die RechengréBen nur durch 
endliche Dezimal- oder Dualbriiche dargestellt werden, und zwar ist die tibliche 
Genauigkeit bei programmgesteuerten Rechenmaschinen zirka 12 Dezimal- 
oder 40 Dualstellen. Das Vorzeichen erfordert bei beiden Systemen eine Dual- 
ziffer, indem beispielsweise + als 0 und — als L geschrieben wird. Was aber 
die Kommastellung anbelangt, so unterscheidet man zwei Haupttypen: 

A. Maschinen mit beweglichem Komma (floating decimal [binary] point), zu 
denen der Bell Computer [6], Mark II [3] und das Rechengerat von ZusE [66] 
gehéren und welche alle Zahlen in sogenannter halblogarithmischer Form 
darstellen: x = + pB*4 (B = Basis des Zahlsystems). Damit ist jede Zahl 
durch ein Zahlenpaar vertreten, z. B. bei Mark II gemaB Fig. 3. 


Se ea ae 
x =[+] 1,7543 85965 | —|2| 
t 7 ial 
1 40 1 


4 = total 46 Dualstellen. 


Fig. 3 
Darstellung einer Zahl in Mark II. 
1) Vgl. BERKELEY [12], S. 133. Im iibrigen ist diese Maschine in [6] beschrieben. 
2) Vgl. [16]. 
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gq heiBt der Exponent von x; die Zahl p, die immer im Intervall 1 = P<B 
liegen muB, wird gelegentlich als Mantisse bezeichnet (in der Tat ist ja ¢ 
die Kennziffer, log p die Mantisse von log” x). 


Sonderwerte. Da die Zahl 0 keine halblogarithmische Darstellung besitzt, — 
muB sie durch ein besonderes Zeichen markiert werden, zum Beispiel sind bei 
Mark II in diesem Fall alle 46 Dualziffern = 0. Der deutsche Konstrukteur 
ZuseE hat auBer fiir 0 auch fiir co sowie fiir die ganzlich unbestimmte GréBe 
0/0 (= co — ow oder dergleichen) je ein besonderes Zeichen eingefiihrt. K. ZUSE, 
der die Zahlen in seiner Maschine durch 32 Dualziffern darstellt’ (24 fiir p, — 
8 fiir g), fiigt zu diesem Zweck noch eine weitere Dualziffer mit folgender Be- — 
deutung bei: Wenn diese 0 ist, so bestimmen die tibrigen 32 Dualziffern eine © 
normale Zahl x = p B4; ist diese aber L, so liegt ein Sonderwert vor, und die ~ 
iibrigen 32 Dualziffern entscheiden, ob es sich um 0, oo oder ? (Unbestimmt) 
handelt. 


B. Maschinen mit festem Komma (fixed decimal [binary] point). Bei diesen 
werden die Zahlen ebenfalls als endliche Ziffernfolgen dargestellt, aber im 
Gegensatz zu den Maschinen mit beweglichem Komma ist hier dem Komma _ 
innerhalb der Zifferngruppe eine feste Position zugewiesen, namlich bei den 
meisten Maschinen am Anfang der Zifferngruppe!) (so daB alle Rechen- 
groBen im Intervall — 1< x< 1 liegen miissen), wahrend sie bei Mark I 
und Mark III wenigstens noch fiir jedes Problem wahlbar ist: 


~ cmt’ ate nay gi 


Mogliche Kommastellungen bei Mark III: ,7,3,6 2,51,4039281706, 
ebenso bei BINAC und IAS: OL L0ULLOOLE00ULLUOLESs- 


Das Vorzeichen erfordert ebenfalls eine Dualziffer, eine Angabe des Expo- 
nenten ist unnotig. 


Fiir den Uneingeweihten scheint diese starre Behandlung des Kommas 
nicht nur ein schwerer Nachteil, sondern iiberhaupt ein unhaltbarer Zustand 
zu sein. Tatsachlich ist es aber so, da8 man, falls die GréBenordnung der in der 
Rechnung auftretenden Zahlen einigermaBen bekannt ist, das Problem so 
transformieren kann, daB dann alle in die Rechnung eingehenden Gré8en im 
Intervall —1< x <1 liegen. Man wird dabei vor allem darauf ausgehen, eine 
nach oben sichere Abschatzung zu erhalten, selbst auf die Gefahr hin, daB 
nachher die in der Rechnung auftretenden Zahlen sehr viel kleiner als 1 sind 
(dafiir hat man ja bei diesen Maschinen zum Teil so viele Stellen). 

AuBerdem hat es der das Problem Vorbereitende als ultima ratio in der Hand, 
die Rechnung so zu gestalten, daB die Skala — welche sich bei den Maschinen 


mit beweglichem Komma automatisch anpaBt — ebenfalls beweglich wird, 
vel. [15]. 


1) Es scheint allerdings, da® die Stellung des Kommas zwischen der ersten und zweiten Ziffer 
vorteilhafter ware, vgl. [35]. 
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Zusammenfassend kann gesagt werden, daB zwar der Umstand, da8 man 
bei einer Maschine mit festem Komma die GréBenordnung der in der Maschine 
gebildeten Rechenresultate ungeféhr vorausbestimmen mu8, um erfolgreich 
damit arbeiten zu kénnen, ein entschiedener Nachteil ist, der aber durch eine 
Anzahl von Vorteilen gegentiber den Maschinen mit beweglichem Komma 
mehr als ausgeglichen wird. Zweifellos wird man in der Weiterentwicklung ver- 
suchen, die Vorteile beider Systeme in einer Maschine zu vereinigen. 


3.4. Rechnen mit halblogarithmisch dargestellten Zahlen 


Es soll hier kurz erklart werden, welche arithmetischen Probleme die vier 
Grundoperationen in einer Maschine mit beweglichem Komma aufwerfen. Be- 
treffs schaltungstechnischer Probleme vergleiche man die sehr ausfiihrliche 
Beschreibung von Mark II [3]. 


Addition und Subtraktion 


Zunachst ist klar, daB nur Zahlen mit gleichen Exponenten addiert werden 
kd6nnen. Deshalb muB man zuerst die Differenz der Exponenten bilden und 
dann den Summanden mit dem kleineren Exponenten um eine gewisse Anzahl 
Stellen nach rechts verschieben, zum Beispiel bei einer dreistellig dezimalen 


Maschine: 
a) 5,32 10-4 0,05 x 10-2 
4,73 10-2—> 4,73 x 10-2 


Korrekte halblogarithmische Summe 4,78 X10-2 


b) 9,82x102 + 9,82x 102 


Die Mantisse ist hier gréBer als B = 10, deshalb hat ein Ausrichten der Summe 
zu erfolgen, deren halblogarithmische Darstellung lautet: 1,06 x 103, 


c) — 5,55X 101 
ea LOA 


— 0,02 x 101 


Auch hier geniigt die Mantisse nicht mehr der Bedingung 1S p< B, es ist 
jetzt eine Stellenverschiebung nach links notwendig, die korrekte Darstellung 
wird so: -2,00x10-!. 


Im Beispiel c erscheinen infolge der Stellenverschiebung nach links am 
Ende eine Anzahl von Nullen, welche eine nicht vorhandene Genauigkeit vor- 
tiuschen. Es ist dies ein Nachteil der Maschinen mit beweglichem Komma, daB 
sie den beim Subtrahieren fast gleich groBer Zahlen auftretenden Genauig- 
keitsverlust nicht erkennen. 


Multiplizieren und Dividieren 


Fiir diese beiden Operationen addiert bzw. subtrahiert man die Exponenten 
und multipliziert bzw. dividiert die Mantissen. 
k 


Nun kann aber das so erhaltene Mantissenprodukt eine Zahl zwischen 1 une 
B? sein, so daB eventuell eine Rechtsverschiebung derselben und gleichzeitig eine ~ 
Erhéhung des Exponenten nétig ist. Entsprechend kann der Mantissenquotient ~ : 
zwischen B-1 und B liegen und verlangt also eventuell nach einer Linksverschie- — 
bung mit Verminderung des Exponenten. Im Detail werden diese beiden Opera-~ 
tionen ahnlich ausgefiihrt wie bei den Maschinen mit festem Komma, fiir welche ' 


sie ausfiihrlich beschrieben sind (vgl. §§ 3.52 und 3.53). : 


? 


Wie man sieht, hat die halblogarithmische Darstellung der Zahlen in einer 
Rechenmaschine den Nachteil, daB Addition und Subtraktion, welches bei — 
weitem die haufigsten Operationen sind, verhaltnismaBig kompliziert werden. 
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Rechnen mit Sonderwerten 


Der Exponent in der halblogarithmischen Darstellung einer Zahl kann nur’ 
in gewissen Grenzen variieren (z. B. von —15 bis +15 bei Mark II), so daB 
die GréBenordnung der in der Maschine vorkommenden RechengréBen auch 
nicht ganz unbeschrankt ist. Es kann nun vorkommen, daB bei der Ausfithrung 
der Rechnung Zahlen entstehen, deren Exponenten die erwahnten Grenzen 
uberschreiten, d. h. der Exponent iiberflieBt?). In einem solchen Fall muB eine 
Rechenmaschine normalerweise anhalten, da sonst falsche Resultate entstehen. 
Das Rechengerat von K. Zuse bildet hingegen einen Sonderwert «sehr groB» ~ 
und rechnet mit diesem weiter, solange dies auf eindeutige Weise méglich ist. 
Erst wenn Operationen wie 0 -co und co — co zu unbestimmten Resultaten 
fithren, kommt es zum Unterbruch der Rechnung. Das unvorhergesehene An- 
halten der Maschine wird dadurch nicht nur verzégert, sondern in vielen Fallen 
uberhaupt verhindert. Zum Beispiel kommt es beim Auswerten von Ketten- 
briichen vor, daB ein Teilnenner co oder sehr groB wird. Da diese GréBe aber 
in einem Nenner steht, ist sie véllig unschddlich; tatsachlich wird auch das 
Rechengerat von ZusSE in einem solchen Fall (im Gegensatz zu andern Ma- 
schinen) nicht anhalten, sondern voriibergehend den Sonderwert co oder «sehr 
groB» bilden. 


3.5. Rechnen mit Maschinen mit festem Komma?) 


Es wird im folgenden immer angenommen, da8 in einer solchen Maschine 
das Komma so gesetzt sei, daB alle Zahlen zwischen — 1 und + 1 (Grenzen aus- 


1) UberflieBen hei®t: Die héchste Stelle des betreffenden Zahlwerks springt von B—1 auf 0 
und der resultierende Ubertrag geht verloren (wie es auch bei Biirorechenmaschinen beobachtet 
werden kann). Es resultiert also ein Fehler von B Einheiten der héchsten Stelle. 

2) Vel. hieriiber auch Burks, GOLDSTINE und vy. NEUMANN [20], § 5, sowie A. W. Burks 
(Lecture 8 in [89]). 
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geschlossen) liegen, d. h. sie treten als Ziffernfolgen fester Stellenzahl N auf, 
bei denen das Komma stets vor der ersten Stelle zu denken ist, zum Beispiel 


B141592003 23 —— = og mit N=10, 


Benn auch bei Mark I und Mark III “ Kommastellung innerhalb oder am 
Ende der Ziffernfolge gewahlt werden kann, so ist dies nur eine Rechenerleich- 
terung und andert das Prinzip nicht. 


3.51. Addition und Subtraktion, Komplementbildung 


Vor Beginn jeder Rechenoperation befinden sich die beiden Operanden 
(Augend und Addend oder Minuend und Subtrahend) im Speicherwerk der 
Maschine. Da sie selbst noch mit Vorzeichen versehen sind, kann erst nach dem 
Ablesen derselben entschieden werden, ob effektiv eine Addition oder Subtrak- 
tion auszufiihren ist. Diese Entscheidung wird dadurch vereinfacht, daB man 
Addition und Subtraktion gleichwertig als Addition von Zahlen behandelt, von 
denen eine oder beide negativ sein kénnen, und auBerdem negative Summanden 
im Rechenwerk nicht als Absolutbetrag mit Vorzeichen, sondern als Komple- 
mente (konegative Gestalt) darstellt: 

Zu jeder negativen Zahl x, die ins Rechenwerk eingeht, wird automatisch 
eine Konstante C addiert, die Summe x = x + C heiBt dann das Komplement 
von x. Wenn diese Konstante C = 2 — 7 ist (wobei 7 = B~* eine Einheit der 
letzten Stelle bedeutet), so ist wegen x > —1, genauer x 2 —1-+ 7 offenbar 
ge 

Somit besteht dann eine eindeutige Unterscheidung zwischen den positiven 

Zahlen und den in komplementarer Form dargestellten negativen Zahlen; 
diese sind unechte, jene echte Dezimal- oder Dualbriiche. Diese Unterscheidung 
ist natiirlich unerlaBlich, so daB nur Konstanten C = 2 — 7 brauchbar sind. 
Da also die im Rechenwerk auftretenden konegativen Zahlen (negative Zahlen 
in komplementiarer Form) = 1 sind, muB dieses noch mit einer zusatzlichen Stelle 
vor dem Komma (Einer) versehen sein. 
' Wird auBerdem fiir positive Zahlen * = x gesetzt, so ist immer * = % 
(mod C), damit auch * + y = x+y (mod C). Um also die richtige Darstellung 
fiir ++ y zu erhalten, miissen ¥ und y addiert und die Summe mod C ins 
Intervall 0 < *%<C oder 0< x SC) reduziert werden. Damit ist das Pro- 
blem der allgemeinen Addition offenbar gelést. 

Von den vielen méglichen Werten fiir C sind eigentlich nur zwei wesentlich 
verschiedene im Gebrauch, namlich B und B — », und man unterscheidet dem- 
gemaB auch zwei Arten von Komplementen, namlich die B- und die (B — 1)- 
Komplemente: 


¥ 


1) Die erste Ungleichung entspricht dem Fall, daB 0 als positive Zahl dargestellt wird, so daB 
0 = 0 ist, wahrend die zweite Ungleichung in Kraft tritt, wenn 0 als negative Zahl behandelt wird, 


wobei dann 0 = C wird. 
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Bei den B-Komplementen (C = B) wird die Null immer als positive Zahl auf- 
gefaBt. Das Rechnen mod B wird dadurch realisiert, daB der Ubertrag von der 
héchsten Stelle (Einer), der ja ins Leere geht, unterschlagen wird; damit bleiben 
alle Zahlen automatisch im IntervallO S$ «<< B=C. 


Beispiel: (Dezimal) : 


+ ,6857—> 0,6857 — ,3456—> 9,0544 

— ,3928—> 9,6072 — ,4321> 9,5679 

fallt weg <— 0,2929 fallt weg <(1)o,2223 
also * + y = 0,2929 Also. yi Lees 
“+y=—-—0,7777 


Bei den (B—1)-Komplementen ist C= B—y, die Null wird dann zweck- 
maBig als negative Zahl betrachtet und folglich durch B— 7 (alle Ziffern sind 
B—1) dargestellt. Das Rechnen mod B— wird so realisiert, daB man den 
von der héchsten Stelle (Einer) ausgehenden Ubertrag, der den Wert B hat, der 
letzten Stelle, deren Einheit den Wert 7 hat, zufiihrt. Dies ist der sogenannte 
Endiibertrag (end around-carry), er entspricht offenbar einer Subtraktion von 
B= 7 

Beispieles (Bit) =e 


+ ,6857—> 0,6857 — ,3456—> 9,6543 

— ,3928—> 9,6071 — ,4321—> 9,5678 
0,2928 1. j932291: 

Enditibertrag \—» 1 Endiibertrag \—» . 1 
Summe 0,2929 Rope | EOD ae 


Die Zuteilung der Null zu den negativen Zahlen bei Maschinen mit (B — 1)- 
Komplementen hat folgenden Grund: 

Eine Summe von positiven Zahlen (auch konegative Zahlen sind ja positiv) 
kann nie 0 sein, so da8 die Null als Rechenresultat im allgemeinen zunachst als 
C und nicht als 0 erscheint. Da aber bei C= B— » der Endiibertrag ausbleibt, 
erfolgt keine Reduktion mod C auf 0, sondern die Null bleibt als 0 = B— 
im Rechenwerk und tragt damit alle Merkmale einer negativen Zahl. 


Komplementbildung 


\ 


Wenn eine negative Zahl aus dem Speicherwerk als Summand ins Rechen- 
werk iibergefiihrt wird, so muB eine Komplementbildung erfolgen, und wenn 
man eine negative Summe speichern will, so muB wieder der Absolutbetrag mit 
Vorzeichen gebildet werden, wozu eine neue Komplementbildung erforderlich 


ist). Es ist deshalb auBerordentlich wichtig, eine einfache Regel fiir die Komple- 
mentbildung zu haben. 


) Einige Maschinen, namlich Mark I [1] und IAS [20], verwenden allerdings die konegative 


Form der negativen Zahlen auch im Speicherwerk, was dann aber Multiplikation und Division 
kompliziert. 
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Wie man sich leicht tiberzeugt, erhalt man das (B — 1)-Komplement einer 

negativen Zahl einfach dadurch, daB man alle Ziffern ihres Absolutbetrages 
zu B — 1 erganzt. Im Dualsystem ist dies besonders einfach: Jede 0 ist durch 
Lund jede L durch 0 zu ersetzen, was auch technisch sehr leicht durchzufiihren 
ist. Falls aber die Verschliisselung der 10 Ziffern des Dezimalsystems geeignet 
gewahlt wird, genieBt man dieselben Vorteile auch in diesem (vgl. Postulat C, 
By 3.2). 

Um hingegen die B-Komplemente zu bilden, ist eine richtige Rechenopera- 
tion erforderlich; man kann namlich zunachst das (B — 1)-Komplement bilden 
und dann in der letzten Stelle eine 1 addieren. Obgleich dies eine sehr einfache 
Addition ist, mu8 doch damit gerechnet werden, daB sie eine ganze Reihe von 
Ubertragen auslist, weshalb die B-Komplementbildung wesentlich mehr Auf- 
wand erfordert. 


3.52. Multiplikation 


Wenn negative Zahlen als Absolutbetrag mit Vorzeichen im Speicherwerk 
sind?), so werden Multiplikator (MP) und Multiplikand (MC) dem Rechenwerk 
in folgender Form zugefiihrt: 


N N 
MPI eS xg Bo* MCS yep)" 9, Bo. 
1 1 


Dabei kénnen die x, und y,, die Ziffern der beiden Faktoren, nur die Werte 
0, 1, ..., B—1 annehmen. Die Vorzeichen der beiden Faktoren werden im 
folgenden immer auBer Betracht gelassen, weil die Maschine nur die Absolut- 
betrage multipliziert und nachher noch das Vorzeichen des Produkts bestimmt. 

Die eigentliche Multiplikation wird nun so ausgefiihrt, daB jede Stelle des 
Multiplikators mit dem ganzen Multiplikanden multipliziert und das Ganze 
summiert wird, d. h. man bildet: 


ie 
y= Sa, y BH®. 
1 


Die durch ein Beispiel (B = 10, N = 4) angedeutete altbekannte Methode: 


pli 79999 = 09999 
9999 


, 11108889 
ist in programmgesteuerten Rechenmaschinen, die mit Relais oder Elektronen- 


1) Wenn die negativen Zahlen auch in komplementarer Form gespeichert werden, so mtissen 
erst — wie bei Mark I — die Absolutbetrage gebildet werden, oder man kann auch direkt mit den 
eventuellen konegativen Zahlen zur Produktbildung schreiten. Hiertiber Naheres in [20], Abschnitt 


5.10 und 5.11, S. 16ff. 


a 
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réhren arbeiten, undkonomisch, weil es ein (2 N)-stelliges Addierwerk braucht, um 
das Produkt zu bilden, wahrend das in dieser Hinsicht sparsamere, ebenfalls wohl- 
bekannte Verfahren der abgekiirzten Multiplikation zu anderen Unzulanglich~ 

keiten fiihrt?). 


Hingegen kann man bei Anwendung des nachstehend beschriebenen Algo- — 
rithmus [Formel (3.5)] mit einem N + 1-stelligen Addierwerk AC?) und einem 
N-stelligen Register MR auskommen (das Register MR muB nicht addieren 
kénnen und ist deshalb konstruktiv wesentlich einfacher als AC) : 


: 
| 


| 
, 
. 
1 
rf 
? 


py = 0 (AC wird vor Beginn der coe geléscht) Bos. 
bra=zy (bet ey) fiir RaSh eh het 
XV=PDo [a 
Realisiert wird dies dadurch, daB man die Stellen B®, B-1, ..., B~* “= 


p, (R= N, ..., 0) im Akkumulator AC und die Stellen B™ N-1, 5 sine eae 


im Register MR aufbewahrt. Die Bildung von #;_, = B-+ (pf, + x, = Ba 


dann, da8 zunachst zu p; ein ganzes Vielfaches von y, welches auch ein ganzes 
Vielfaches von 7 = B~* ist, addiert wird. Dieser Vorgang spielt sich allein in 
AC ab und bewirkt keinerlei Anderungen in MR. Alsdann hat noch eine Multi- 
plikation mit B-1, das heiBt eine Verschiebung um eine Stelle nach rechts, und 
zwar gemeinsam in AC und MR, zu erfolgen; die letzte Stelle von AC wird 
dabei zur héchsten in MR. 


Wir zeigen diesen Vorgang am Beispiel * x y = ,2345 X ,6789: 

AC MR (MR) 
Pa 0.0000 He (2345) 
Prt 5y=10f; 33945 eee (2345) 
Ps 03394 5x (5234) 
pa t4y=10, 30550 5x (5234) 
Ps 03055 O5** (0523) 
Pet 3y=10, 23422 O5** (0523) 
py 02342 205* (2052) 
Pit 2y=10 py 15920 205* (2052) 
pPy= x= 01592 0205 (0205) 


Wie man sieht, sind die mit einem Stern markierten Stellen in MR fiir das 
Ergebnis belanglos, da sie infolge der fortlaufenden Rechtsverschiebungen all- 
mahlich eliminiert werden. Man kann deshalb am Anfang der Multiplikation 
die N Stellen des Multiplikators in MR speichern und hat dann die fiir die Be- 
rechnung von #;_, bendtigte Ziffer x, des Multiplikators immer genau in der 
letzten Stelle von MR zur Verfiigung. Was das Register unter diesen Umstanden 
enthalt, ist in der obigen Tabelle rechts auSen in Klammern angegeben. 


1) Das Rechnen mit gréGerer Genauigkeit (§ 3.55) ist unméglich. 
2) Dieses wird natiirlich auch fiir die in § 3.51 behandelte Addition verwendet. 
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_ Nach ausgefiihrter Multiplikation steht in den beiden Registern AC und 
MR das (2 N)-stellige Produkt der beiden N-stelligen Faktoren zur Verfiigung, 
namlich die ersten N (sogenannte wesentliche) Stellen in AC, die letzten N 
(sogenannte wnwesentliche) Stellen in MR. Normalerweise werden nur die 
wesentlichen Stellen des Produkts verwendet (dieses wird also auf N Stellen ab- 
gerundet, vgl. § 3.6), denn samtliche N unwesentlichen Stellen sind als ungenau 
zu betrachten, wenn die beiden Faktoren mit Fehlern der GréBenordnung 7 
behaftet sind; sie werden jedoch fiir das Rechnen mit hdherer Genauigkeit 
(§ 3.55) aktuell. 

Die besprochene Methode [Formel (3.5)] ist prinzipiell fiir alle Zahlsysteme 
brauchbar. Aber beim Dualsystem wird sie besonders einfach, was einer der 
Vorziige desselben ist. Die Operation 


i 
= (but %e VY) = Pra 
in (3.5) reduziert sich dann auf?) 


{ Bs fiir “x, = 0 | 
204 =" es Bay tg (3.6) 
|p. + 9 fir vaca, | 


Beim Dezimalsystem (und auch bei allen andern Zahlsystemen mit B a) 
ist dagegen x, y im allgemeinen eine wenn auch einfache Multiplikation, die 
auf verschiedene Arten ausgefiihrt werden kann: 

a) Durch fortlaufende Aufsummierung von y wie bei den Biirorechenmaschi- 
nen (Bell-Computer). Man kann die Arbeit dadurch etwas abkiirzen, daB 
man fiir x, =5 von 10 an abwarts zahlt; dies wird insbesondere beim 

Bell-Computer so gemacht. 

Man bildet zuerst nach a die Vielfachen 2y, 3y, 4y, ... bis 9y und spei- 

chert diese in geeigneten Registern. Dann erst beginnt die Multiplikation 

gemaB (3.5), wobei diese Vielfachen verwendet werden (Mark I). 


c) Man bildet nur die Vielfachen 2y, 4¥, 5y, aus denen man zusammen mit 
y in allen Fallen x,y durch héchstens eine Addition erhalt (Mark II). 
d) Durch die Methode der Links- und Rechtskomponenten?) (SSEC, Mark II 
und Lochkartenmaschinen) : 
N 


= 


Die Zahl x, y =D) %, y; B-’ wird in zwei Summanden zerlegt, indem 
1 
man jedes der Elementarprodukte ~;, yj wie folgt aufspaltet: 
Kp Vj= Buy + n;, mit 0 =i) = BB, 


1) Vgl. insbesondere Burks, GOLpDsTINE und NEUMANN [20]. 
2) Vgl. Manual of Mark III [4]. 
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Dadurch wird nun: 


yy = BS tap BI + 3) ty BY = B L(x, y) + R(x, y) 


(zum Beispiel ist fiir x, = 7; y= 7385: L = ,4253; R= ,9165}. Die Links- 
und Rechtskomponenten aller Elementarprodukte 
[0 x 0 bis (B — 1) x (B—1)] 


miissen in einer solchen Maschine natiirlich fiir immer gespeichert werden’). 
Gegeniiberstellung der vier Methoden zur Bildung von x, y: Die Anzahl der 
Additionen zur Ausfiihrung der Multiplikation N-stelliger Zahlen betragt im~ 
Dezimalsystem im Mitel: 


im Falle a: 44.N 


a 


1S ea INE 
Se iis 
ole ZN 


Ferner betragt die entsprechende Vergleichszahl im Dualsystem [Multiplikation 
gemaB Formel (3.6)]: 1,66 N, weil im Dualsystem zur Erreichung derselben — 
Genauigkeit gem4B § 3.1 zirka 3,32mal so viele Stellen notwendig sind wie im 
Dezimalsystem. 

3.53. Division und Quadratwurzel 


Wahrend Addition, Subtraktion und Multiplikation Standardoperationen 
sind, die in keiner verniinftigen Rechenmaschine fehlen diirfen, kGnnen nicht 
alle programmgesteuerten Rechenmaschinen dividieren; zum Beispiel bilden 
Mark II, EDSAC und Mark III den Quotienten nicht auf die tibliche Weise, 
sondern durch Multiplikation des Zahlers mit dem auf andere Weise erhaltenen 
reziproken Wert des Nenners?). Die Quadratwurzel wird sogar nur von wenigen 
Maschinen (ENIAC, Bell Computer, Rechengeraét von ZUSE) nach dem bekann- 
ten elementaren Verfahren gezogen, die tibrigen berechnen sie durch Auflésen 
der Gleichung x? = a nach dem Newtonschen Verfahren (§ 4.8). 

Division. Sei vy der Absolutbetrag des Nenners, z derjenige des Zahlers. Da 
nun die Maschine nur mit Zahlen rechnen kann, deren Absolutbetrag kleiner 
als 1 ist, muB die Nebenbedingung y > z erfiillt sein. 

Die zweckmaBigste Art zu dividieren ist die Umkehrung der in § 3.52 dar- 
gelegten Methode der Multiplikation: 


eae Se 
Pr = Bopia—%y (R=1,..,N), (3.7) 


py = BN-facher Rest der Division. 


1) Vel. hierijber § 5.14. 
*) Vgl. Description of a Relay Calculator [3], Manual of Mark III [4] sowie § 4.8 dieser Arbeit. 


* 
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Dabei soll x, die gréBte ganze Zahl mit der Eigenschaft sein, daB noch £, = 0 

ist, dann ist 0,%4%.%3 ++: ¥y der ohne Riicksicht auf die noch folgenden Stellen 
auf N Stellen abgerundete Absolutbetrag des Quotienten. 
_ Realisiert wird der Algorithmus (3.7) dadurch, daB man zunichst p» nach 
AC bringt, eine Stelle nach links schiebt und dann so oft y-subtrahiert, bis in 
AC eine negative Zahl erscheint. Wenn dies geschieht, wird einmal y addiert, 
um , positiv zu machen (Wiederherstellung des positiven Restes). Dann ver- 
schiebt man eine Stelle nach links und subtrahiert wieder y usw. 

Im folgenden Beispiel ist die dreistellige Division 0,352:0,546 ausgefuhrt: 


AC MC (MC) 

a Po 0352 000 (000) 
10 Po fe ie 3520 000 (000) 

Pi i 0244 000 (006) 

10 p, G5 2440 000 (060) 

Pe z 0256 000 (064) 

10 py CUETS, 2560 000 (640) 

“Dy pad 0376 000 (644) 


Also ist 0,644 der abgerundete Quotient und 0,000376 der Rest. Da das 
‘Register MR hierbei unbeniitzt bleibt, kann man es zur Aufnahme der fort- 
laufend gebildeten Ziffern ~,, des Quotienten beniitzen, was rechts nebenste- 

hend in Klammern ausgefithrt ist (vgl. den entsprechenden Kunstgriff bei der 
~Multiplikation, § 3.52). 

Dieses Verfahren bewahrt sich gut bei Maschinen, welche negative Zahlen 

als B-Komplemente behandeln, aber mit (B — 1)-Ko mplementen ist es besser, 

‘die Formeln (3.7) dahin abzuandern, daB man y und z durch —y und —z 

ersetzt, so da dann «; die gréBte Zahl mit der Eigenschatt ist, daB p, <0 ist?). 

Im iibrigen kann die Methode natiirlich fiir jedes Zahlsystem Verwendung 
finden, aber im Dualsystem 1aBt sich eine erhebliche Vereinfachung erzielen: 


Division ohne Riickstellung des Restes im Dualsystem [20], [67] 
(Non restoring division) 
Der Algorithmus 


Po=2; | 
by = 2 Pera — Y Sgn (Px-1) 
2 xy = 1+ sen (Pr) | 


liefert den gesuchten Quotienten z:y als 0,%1%-.. %y im Dualsystem, wobei 
die Reste p, im Gegensatz zu (3.7) beliebiges Vorzeichen haben kénnen. 


Es gilt namlich, wie man durch vollstandige Induktion zeigen kann: 


pees 
pei (2 ty 2-4) y| (1%) Ys 
1 


(3.8) 
Sr ee | 


1) Vgl. A. SPEISER, Entwurf eines elektrischen Rechengerates, Dissertation ETH. (Ziirich 1950). 


woraus weiter folgt: 

N 3 ee 
1 (Sar*)y =2-"|py+ (1—4y)y|S2-%, 
denn entweder ist Py = 0, damit nach (3.8) #y =1, oder dann haben die beiden 
Glieder py und (1 — xy) y ungleiches Vorzeichen. 


Die Division ohne Riickstellung des Restes kann im Prinzip auch auf das 
Dezimalsystem tibertragen werden, bringt dort aber keine wesentlichen Vorteile. 

Quadratwurzel. Da die bekannte elementare Methode des Quadratwurzel- — 
ziehens nur eine Division mit variablem Divisor ist, kann man sie sehr wohl auf — 
programmgesteuerte Rechenautomaten iibertragen, Ein wirklich schénes Ver- 
fahren erhalt man aber nur im Dualsystem, wenn man auf die Riickstellung 
des Restes verzichtet!), der Algorithmus, welcher die verzifferte Quadrat- — 
wurzel 0,%,%» ... Xy von z liefert, lautet dann: 


Pe = Cees 
/ k-1 : 
pr=2 Pas = sgn @xa)(2*+ 37 xj 24|+ ‘a i he 


2 x, = 1+ sgn (p;) - 


Beweis: Wir zeigen vorerst, daB fiir k > 0 
” Rd \2 
py = 2 E (2* + D4; 2-4) | ' (3.10) 
1 ) 


was zunachst ftir k = 1 richtig ist, da wegen z > 0 immer Pp, 
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= 2z—1/2 ist. 

k-2 2 

Ferner folgt mit der Annahme p,_, = 2"? E _ (2- ae Sy je) | aus (3.9): 
1 


k-1 


In beiden Fallen: pz., = 0(%,_,=L) und p, ,<0(%,_,= 0) kann man diesen. 
Ausdruck vereinfachen und zeigen, daB er mit (3.10) tibereinstimmt. Damit 
ist (3.10) richtig, und man kann daraus ablesen, da8B 


N 2 
z— (z Xn zt) 
1 


3.54. Skalarfaktoren 


Beim Rechnen mit einer Maschine mit festem Komma miissen alle in der 
Maschine auftretenden Zahlen im Intervall —1< «<1 gehalten werden, was 
nicht leicht ist, da ja Addition und Division aus diesem Bereich herausfiihren. 
Das Uberschreiten oder Erreichen dieser Grenzen beim Addieren nennt man 
das UberflieBen des Addierwerks AC. Zwar hat letzteres eine Einerstelle, so 
daB die Summe in AC noch korrekt gebildet wird, aber beim Speichern ent- 
stehen Fehler, weil die Speicherzellen diese Einerstelle nicht haben. Obwohl die 


= Se7, » Ww.z.b.w. 


1) Diese Methode stammt von K. Zuser [67]. 


7 
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meisten Maschinen in einem solchen Falle Alarm geben, mu8 das UberflieBen 
des Addierwerks unter allen Umstanden verhindert werden (durch sorgfaltige 
Vorbereitung des Problems), weil die Untersuchung der Situation und die Be- 
hebung des Ubelstandes sehr zeitraubend sind: 

Wie schon erwahnt, mu8 ein numerisches Problem eventuell zuerst einer 
Transformation unterworfen werden, damit alle RechengréBen — die zunachst 
noch ganz unbekannt sind — ins Intervall -1<x*<1 zu liegen kommen, 
das heiBt, man gibt eine Zahl x als x* = B-™ x in die Maschine ein, wobei m 
eine durch Abschatzung erhaltene ganze Zahl mit der Eigenschaft | «|< B™ 
ist!). Es empfiehlt sich aber als SicherheitsmaBnahme, den Exponenten m etwas 
-gréBer zu wahlen, als unbedingt notwendig erscheint. 


Die in einem Problem auftretenden Zahlen kénnen nun in Klassen von solchen, 
die mit dem gleichen Skalarfaktor B-™ ins Intervall (— 1, 1) reduziert worden 
sind, eingeteilt werden. Es ist dabei in den meisten Fallen méglich, die Klassen- 
einteilung durch die den Rechengré8en in der physikalischen Formulierung des 
Problems anhaftenden MaBeinheiten und Dimensionen vorzunehmen. Es mtissen 
dann kaum jemals GréBen verschiedener Klassen addiert werden, so daB die 
Addition keine Probleme aufwirft (andernfalls hatte man mit ahnlichen Kompli- 
kationen zu rechnen wie bei Maschinen mit beweglichem Komma). 

Bei Multiplikation und Division sind jedoch Stellenverschiebungen die 
Regel: 

Es sollen z= xy und w= x/y gebildet werden, wobei die vier GréBen 
Pw, 2, v bzw. alsa* = BB? x, y¥ = B-*y, z* = B-" z und w* = B- w in der 
Maschine seien. Dann ist: 


z* = Br xy = B- B? Ba(x* y*) = BPte-r(x* y*), 
ney ib ad Ek Same ot a : 
w= B (=) = B-sB?B dberd == pre (=). also: 


Nach jeder Multiplikation ist noch eine Verschiebung des Produkts um 
p+4q-—v Stellen nach links notwendig?), ebenso nach jeder Division eine Ver- 
schiebung um ~ — g — s Stellen nach links. 

Dabei sind 4, g, 7, s die Exponenten der Klassen, denen x, y, z, @ angehoren. 


3.55. Rechnen mit hiherer Genauigkert*) 


Die programmgesteuerten Rechenmaschinen, soweit sie mit festem Komma 
arbeiten, kénnen trotz ihrer begrenzten Stellenzahl mit beliebig groBer Genauig- 
keit rechnen, sofern man den Mehraufwand in Kauf zu nehmen gewillt ist. Als 
extremes Beispiel sei erwahnt, daB die Zahlen e und x auf der zehnstelligen 


_ 3) Zahlen, die ihrer Natur nach ganz sein miissen, versieht man zweckmaBig mit dem Skalar- 
faktor 7 = B-N; dies lauft darauf hinaus, da8 man sich in diesen speziellen Fallen das Komma 
in der Maschine hinter der letzten Ziffer denkt. 

2) Unter einer Linksverschiebung um — ” Stellen ist natiirlich eine Rechtsverschiebung um 
n Stellen zu verstehen.. AuBerdem sind Stellenverschiebungen immer von der in § 3.52 beschrie- 
benen Art, d. h. gleichzeitig in AC und MR. 

8) Vgl. [15] und [24], Bd. 1, § 9.8ff. 
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ENIAC in zirka 80 Stunden auf je 2000 Dezimalen genau berechnet worden 
sind [46]. : 
Das Rechnen mit (mN)-stelligen Zahlen auf einer N-stelligen Maschine 


verlauft nun nach folgendem Prinzip: Die zur Basis B verschliisselte Zahl 
nN ; 
x = >| x, B-® wird vom Komma aus in Gruppen von je N Ziffern, das heiSt in 


1 
(nN)-stellige Zahlen &; aufgespalten, was einer Aufspaltung der obigen Summe 
nach aN Eke ; 


vm SB i= = SEB) 


entspricht, wobei die €; ganze Zahlen mit 0 = Ee B® sind. Zum Beispiel — 
Ni tO ,00 = 2 30, 73529: Shee AS, ah ZONE 
Man kann die €; auch als «Ziffern» der Zahl x in einem Zahlsystem mit der ~ 
Basis B™ auffassen. In der Maschine wird nun eine solche Zahl x durch ihre 
(nN)-stelligen «Komponenten» &,, &,...,&, dargestellt und bendtigt deshalb © 
n Zellen zum Speichern. Das Rechnen mit solchen Zahlen x, y, ... muB auf das 
Rechnen mit ihren N-stelligen Komponenten &,;, 7;, ..., welche die Maschine © 
verarbeiten kann, zuriickgeftihrt werden und soll im folgenden behandelt wer- | 
den. Es ist allerdings zu beachten, daB die Komponenten durchwegs ganze 
Zahlen sind, so daB in diesem Abschnitt das Komma bei allen Zahlen und auch 


in den Zahlwerken immer am SchluB der Ziffernfolge ausgefaBt werden soll. 


Addition: z= x*+y= 2) (&;+;) BNi= Se B-Ni, Weil nun é;+ 7; = BY 
1 af 


sein kann, ist nicht einfach ¢; = €; + 7;, sondern es sind gem&8 nachstehendem 
Algorithmus noch Gruppeniibertrage v; zu beriicksichtigen: 


=: (? | 

C; = €; + 95 +; (mod B”) ey, ie ku! (3.11) 
Vj, = ganzer Teil von B~Y (£; + y;+,;) | WS? Baar obeyed | 
vj, ist der bei der Addition von &; + 7; + v; 1m Rechenwerk entstehende 
Uber fluB?). 


Die Subtraktion wird natiirlich durch Komplementbildung auf die Addition 
zuriickgefiihrt, zum Beispiel (N = 2, 1 = 3): 


+ ,529373—> 0’52 0’93 0°73 — ,529373> 9°47 0’06 0’26 
— ,352589- 9’ 64 0’ i4 0710 + ,352589-—> 0’35 0’25 0’89 
om 16 re: 67 0’83 9°82 OFS (as 
se 
il 
~ FEA : 
+ ,176784 <0'°17 0’67 0’84 — ,176784 < 9’82 O32 0’15 


1) In diesem speziellen Fall sind Uberfliisse toleriert und diirfen die Maschine nicht zum Anhalten 


bringen, nur diirfen sie nicht verlorengehen, sondern miissen zur nachst hdheren Gruppe addiert 
werden. 
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(Die Stelle, wo das Komma im Rechenwerk normalerweise zu denken ist, ist 
durch ein hochgestelltes Komma angedeutet). Wie man sieht, werden nur in der 
vordersten Komponente normale (B — 1)-Komplemente (mit dieser Kommastel- 
lung Ergaénzung auf BN+1— 1) verwendet, in den hinteren Gruppen wird auf 
BN — 1 erganzt. 


Multiplikation). Dasexakte Produkt der beiden positiven Zahlen x= Dy oo 4 
und y= ay n; B-? ist eine (2N)-stellige Zahl, von der aber wie bei der ge- 
‘wohnlichen Multiplikation meist nur der aus den ersten a Stellen bestehende 
eeil ae ¢, B-*® gebraucht wird. Es ist nun xy = is 5 E,; n, BONG, wobei 


‘man ais den «Elementarprodukten» €; 4, ganz leg vorgeht wie bei der 


Methode der Links- und Rechtskomponenten, indem man sie wie folgt zerlegt: 
E; My = BY Aj, + Oix- 


Dabei sind die ganzen Zahlen /,;;, und o;;, der wesentliche bzw. unwesentliche 
Teil des Produktes &; 7;,. Man erhalt so schlieBlich die Darstellune 


Vo —— 0) ; 
j 25 i recoe As 2 Q;,,-; + % (mod B”) (3.12) 


_ = ganzer Teil von BX ()" A, 5, + 3) Gs-5+ v;); 
fiir die (nN)-stellige Zahl z* = 3) ¢, B-™', die, wie man zeigen kann, bis auf 
i 


einen Fehler von hichstens 2 1 — 1 Einheiten der letzten Stelle mit dem Pro- 
dukt z = xy iibereinstimmt. Der Algorithmus (3.12) vereinfacht sich im weitaus 
meistgebrauchten Fall der doppelten Genauigkeit zu 


Pr=AotaAnit mod BY 
Se 12 21 O11 ( ) (3.13) 


—_—— 


£,=4n + ganzer Teil von B™* N (Ais Sida + On) - 


Beispiel (N = 2, » = 2): ,7175 x 4492 = ,32230100 (exakt). 
Ase92 = 6532, ° Ag = 65 
74X44 = 3300, Ag, = 33 


Thx 44 = 13124, bon = 247 Ay! 
rape 
eS ye alg 
f, = 22 €, = 32 also 2* = ,3222 


Die Division mit mehrfacher Genauigkeit wird durch Multiplikation des 
Zahlers mit dem reziproken Wert des Nenners bewerkstelligt, welch letzterer 


=) Vel. auch W. EBERL, S. 411 dieses Heftes. 
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durch ein Iterationsverfahren bestimmt wird (vgl. § 4.8), wobei der reziproke 
Wert der héchsten Komponente als erste Approximation dient. Entsprechendes — 
gilt fiir die Quadratwurzel. 

Es ist ganz klar, daB das Rechnen stark verlangsamt wird, sobald die 
Stellenzahl der Maschine nicht mehr geniigt und man deshalb mit mehrfacher 
Genauigkeit arbeiten muB. v. NEUMANN [24] Bd. 1, Kap. 9, zum Beispiel rechnet 
mit achtfacher Rechenzeit fiir doppelte Genauigkeit. 


3.6 Aufrunden 


Die programmgesteuerten Rechenmaschinen réchnen mit Dezimal- oder | 
Dualbriichen fester Stellenzahl N. Da nun gewisse Rechenoperationen — ins- — 
besondere Multiplikation und Division, bei beweglichem Komma auch die Ad- © 
dition — aus diesem Zahlenbereich herausfiihren, mu8 man die tiberschtissigen © 
Stellen weglassen, womit die Resultate verfalscht werden. Es sind folgende © 
Aufrundungsregeln im Gebrauch (vgl. auch Fig. 4)?): 


Maximaler Mittlerer Mittleres 
Fehler Fehler Fehlerquadrat — 

| | | = 

a Feet Lg hg | 7/2 = 
| | ) 

| | | 

. | | | 7? 
? 
~—<—<— ites cos operons [pte y/ 2 0 12 


~ 


} 19 


Fig. 4 
Aufrundungsregeln. 


a) Die tiberschtissigen Stellen fallen einfach weg. 


b) Die letzte nicht weggelassene Stelle wird um 1 erhdht, falls der weggelassene 
Teil 2 7/2?) ist, d. h. man addiert zuerst 1/2 und handelt dann nach a. 


c) Die letzte nicht weggelassene Stelle wird gleich der nachstliegenden unge- 
raden Ziffer gemacht, es sei denn, der weggelassene Teil sei genau 0. Dies 
heiBt sowohl bei einer rein dualen als auch bei einer dezimalen Maschine, 
welche die ungeraden Ziffern durch ungerade Tetraden verschliisselt, daB 


die letzte Dualziffer der Zahl L zu machen sei, ohne die tibrigen zu ver- 
andern. 


1) Vel. Abschnitt 5.12 bei [40]. 
*) 1 = eine Einheit der letzten Stelle (wie in § 3.51). 
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Es ist nun sehr wichtig, daB die Rundungsregel symmetrisch sei, d. h. daB 
der mittlere Fehler verschwinde, weil dieser sich besonders gefahrlich auswirkt. 

Infolgedessen scheidet die Rundungsregel a aus. Da unter den symmetri- 
schen Regeln c vor b den Vorteil hat, daS man den weggelassenen Teil nicht 
kennen muB und daB das Aufrunden keine Ubertrage auslést, wird sie vor allem 
beim Dividieren verwendet (bei Mark II fiir alle Operationen). 

Fiir konegative Zahlen — insbesondere mit (B — 1)-Komplementen — wird 
Aufrunden komplizierter!), so da8 man meist nur die Absolutbetrage rundet. 

Beispiel fiir Regel b (Runden auf zwei Stellen nach dem Komma): 


Zehnkomplement | Neunkomplement 
= 20353-97765 | * 9,7764 
kn > +50 | —50 
—0,22 <9,78(15) | 9,77(14) 


Wahrend man also bei B-Komplementen immer /2 vor dem Weglassen ad- 
dieren mu8, muB man im Fall der (B — 1)-Komplemente bei negativen Zahlen 
7/2 subtrahieren. 


3.7 Ubersetzen vom Dezimal- ins Dualsystem und umgekehrt 


Da das Dualsystem bei den rein dualen Maschinen (§ 3.1) eine rein interne 
Angelegenheit ist, miissen die von auBen kommenden und deshalb im Dezimal- 
system gegebenen Anfangswerte eines Problems zuerst ins Dualsystem umge- 
rechnet werden; auBerdem hat am Ende der Rechnung eine Riickiibersetzung 
ins Dezimalsystem zu erfolgen. 

Um hierzu nicht noch ein spezielles Rechengerat bauen zu missen, 1aBt 
man diese Ubersetzungsarbeit von der Rechenmaschine selbst ausfiihren, wobei 
sich aber eine spezielle Schwierigkeit ergibt: Die Maschine rechnet rein dual, 
muB aber anderseits mit dezimalen Zahlen umgehen kénnen, um sie ins Dual- 
system zu iibersetzen. Diese Schwierigkeit wird dadurch behoben, da man die de- 
zimalen Zahlen verschliisselt, und zwar durch die direkte Verschliisselung (§ 3.2). 
Der Einfachheit halber wollen wir uns im folgenden auf Maschinen mit festem 
Komma beschranken, so daB nur echte Dezimalbriiche zu iibersetzen sind. Sei 


M 
also die Zahl x = . x, 10-* M-stellig dezimal (dual verschliisselt) in der Ma- 
1 


schine. Jede der Ziffern x, ist als vierstellige Dualziffer dargestellt, so daB die 
ganze Zahl 4 M Dualziffern braucht, welches gerade die Kapazitat der Maschine 


sein soll. Wenn man diese 4 M Dualziffern als eine reine Dualzahl auffaBt, 
M 


stellt sie die von x vollig verschiedene Zahl y= D/ x, 16-* dar. Tatsachlich 
: 1 


1) A. Burks, Digital Machine functions, Lecture 8 in [39]. 


— 


~ 
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wird dies von der Maschine auch so aufgefaBt, so daB die Auigabe, vom Stand- 
punkt der Maschine aus gesehen, darin besteht, im Dualsystem aus der Zahl y 
die Zahl x zu berechnen, was durch folgenden Algorithmus geleistet wird?): 


po=Vr W=9, | 
x, = ganzer Teil 
ae : von 16 p;_4”) | 
by, = gebrochener Teil | hz hs Zanes, (3.14). 
Gu 10 Gp lO | 
Rote 
A 0,625¢@ 


Die Zahl 0,625™, welche eine Maschinenkonstante ist, wird natiirlich ein 
fiir allemal in der Maschine gespeichert. 

Die Riickitbersetzung ins Dezimalsystem ist noch wesentlich einfacher*): 
Sei die Zahl x, | « | < 1 im Dualsystem gegeben, gesucht ist dieselbe Zahl dual 
verschliisselt im Dezimalsystem. 


M 
% = DS" 4,10-* (x, = Tetraden) . 
1 


M 
Die Maschine berechnet natiirlich nicht x, sondern wieder die Zahl y= Y) x, 16-*, 
uf 


welche dual geschrieben gleich lautet wie die dual verschltisselte Dezimalzahl x: 
] 


Po mee qo = Ue | | 
X%_, = ganzer Teil 
von 10 £;_4 
py = gebrochener Teil _| lal ees (3.15) 
qu= 16" x, +16 9,4, | | 
Y= Iu | 


Fortsetzung und SchluB8 folgt in ZAMP II, 1951, 


1) Im wesentlichen bei v. NEUMANN und GoLpsTINE [24], Bd. 1, Kap. 9. Eine davon abwei- 
chende Methode wurde fiir die EDVAC gegeben [34]. 

2) Der ganze Teil von 16 Pp Wird praktisch folgendermafen bestimmt: 

Man bildet das Produkt 161~™“ Py—1, dann erhalt man in MR den gebrochenen, in AC 
dagegen den ganzen Teil von 16 Py-1» Wobei in AC das Komma am Ende der Zifferngruppe zu 
denken ist. Denkt man hingegen in AC das Komma normal am Anfang, so ist die in AC stehende 
Zahl gerade die bei der Berechnung von gq, auftretende GréBe 16-M Xp, « 

8) Im wesentlichen bei v. NEUMANN und GoLpstTINE (24], Bd. 1, Kap. 9. Eine davon abwei- 
chende Methode wurde fiir die ED VAC gegeben [84]. 
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Uber die Tragweite der Begriffe «Brennpunkte» und «Brenn- 
‘weite» in der Elektronenoptik und die starken Elektronen- 
linsen mit Newtonscher Abbildungsgleichung 

1. Mitteilung 


Von WALTER GLASER und Otto BERGMANN, Wien?) 


1. Existenz von Brennebenen und Kritik des Brennweitenbegriffes 


Unter dem Brennpunkt eines elektrisch-magnetischen Abbildungsfeldes ver- 
steht man den Achsenschnittpunkt einer achsenparallel einfallenden Elektro- 
nenbahn. Wenn wir iiber den Feldabfall keine Voraussetzungen machen, so 
wird die Lage des Brennpunktes davon abhangen, wo wir die Elektronenbahn 
achsenparallel annehmen. Man denke zum Beispiel an ein homogenes Magnet- 
feld, das sich seiner Natur nach auf beiden Seiten ins Unendliche erstrecken 
muB und in dem daher die Elektronenbahnen bis ins Unendliche ihren oszilla- 
torischen Charakter beibehalten. Infolge der Rotationssymmetrie ist die opti- 
sche Achse stets eine Feldlinie des elektrischen oder magnetischen Feldes. 
Da den Elektronen gerade langs der Achse der Weg freigegeben werden muB, 
kénnen weder Elektroden- noch Polschuhflachen die Achse schneiden. Die 
axiale Feldlinie kann daher weder Anfang noch Ende haben, d.h. die elektro- 
nenoptischen Abbildungsfelder miissen sich grundsatzlich immer nach beiden 
Seiten ins Unendliche erstrecken. Ein abbrechendes Abbildungsfeld, mit dem 
verschiedentlich in der Elektronenoptik argumentiert wird, ist daher eine 
unerlaubte Fiktion. Dies trifft — nebenbei bemerkt — auch fiir ein nur von 
der z-Achse abhangiges Abbildungsfeld zu. Denn das einzige Feld ohne Abhan- 
gigkeit von der radialen Komponente ist das homogene Feld. Die Abbildungs- 
eigenschaften sind aber, zumal im elektrischen Fall, gerade durch die radialen 
- Feldkomponenten gegeben. Rechnerisch ergibt sich aus div€=0 mit 
© = —grad g fiir das Potential p die bekannte Entwicklung 

“ it 
p= 00) — 7 "+. -- 
Ist @” nicht identisch Null, also das Feld nicht homogen, dann muB8 es eine 
radiale Feldkomponente &,= 7/2 @"(z) ++ besitzen. Ein allein von z ab- 
hangiges rotationssymmetrisches Feld kann also nur ein homogenes sein, das 


1) Institut fur angewandte Physik an der Technischen Hochschule Wien und Abteilung fur 
Elektronenoptik der Siemens & Halske AG. in Berlin. 
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sich nach beiden Seiten ins Unendliche erstreckt. Fiir dieses ist namlich 
@" (z) = 0, die Feldstarke kann sich also mit z nicht andern. Hatte es Anfang 
oder Ende, so miiBte an diesen Stellen O"(z) + O sein, was notwendig zu einer 
nichtverschwindenden Radialkomponente fiihren miiBte. 

Ganz analog zeigt man den gleichen Sachverhalt fiir ein magnetisches Feld. 
Die Feldstarken zeigen stets eine Abhangigkeit von der z-Koordinate und haben ~ 
daher nach obigem auch radiale Komponenten. 

Zusammenfassend kénnen wir also sagen, daB (quellen- und wirbelfreie) — 
elektrisch-magnetische Abbildungsfelder stets von der z-Koordinate abhangen ~ 
und sich nach beiden Seiten ins Unendliche erstrecken miissen. 

Da der Brennpunkt der Achsenschnittpunkt der achsenparallel einfallenden 
Elektronenbahn ist, miissen wir uns daher fragen, wie stark das Abbildungsfeld ~ 
nach beiden Seiten abfallen mu8, damit man in entsprechender Entfernung die — 
Elektronenbahn als geradlinig auffassen kann. Die Differentialgleichung der — 
flachen, achsennahen Elektronenbahnen muB daher Lésungen besitzen, die im — 
Unendlichen asymptotisch in Gerade iibergehen. Diese Differentialgleichung 
lautet bekanntlich 


Jb 2 VOo) + (2 + gig Bla) eno. @ 


In (1) bedeutet @ das elektrische Potential langs der. z-Achse, B,(z) ist die 

z-Komponente der magnetischen Induktion auf der Achse. Aus dem linearen 

und homogenen Charakter der Differentialgleichung folgt die Existenz der Ab- 

bildung. Sind @(z) und o(z) zwei unabhangige Lésungen von (1), so ist Dingort 2 

und Bildort z, durch die Abbildungsgleichung?) 
o(Zo) as" G(24) (2) 
Q (Zo) 0(4) 

verkniipft. Die VergréBerung ist durch 


gegeben. Lediglich die Tatsache der optischen Abbildung ist fiir die Anwen- 
dungen von Wichtigkeit, wahrend die engere Analogie mit der geometrischen 
Lichtoptik, wie der Gebrauch der Begriffe Brennpunkt und Brennweite, zwar 
praktisch, aber unwesentlich ist. Im allgemeinen Fall der elektronenoptischen 
Abbildung von diesen Begriffen zu sprechen, wie dies meist geschieht, entbehrt 
uberhaupt jeder Berechtigung. 

Gibt es nur achsenparallel geradlinige Asymptoten, so kénnen wir allein von 
Brennpunkten sprechen. Gibt es auch zur Achse geneigte geradlinige Asym- 
ptoten, dann schneiden sich die mit gleicher Richtung in der durch den Brenn- 
punkt gehenden Ebene, und wir kénnen auch von Brennebenen sprechen. 


1) Vel. zum Beispiel W. Giaser, Z. Phys. 117, 314 (1941). 
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Die Differentialgleichung (1) kann man durch die Substitution 


‘a. 0 = @-i4 R (4) 
auf die Gestalt 

i R' + Q(z)R=0 (5) 
bringen, wobei 


00) = <5 (a) + eae (6) 


ist. Die relativistische Massenkorrektur kann leicht beriicksichtigt werden’). 
Wie im Anhang I gezeigt wird, ist fiir die Existenz des dingseitigen Brenn- 
_punktes die Konvergenz des Integrals 


oo 


fe Q(z) dz (7) 


erforderlich. Fiir die Existenz der dingseitigen Brennebene muB die scharfere 
Bedingung erfiillt sein, daB das Integral 


| # Ql) 4 (8) 
konvergiert. 

Fiir die Existenz des bildseitigen Brennpunktes bzw, der bildseitigen Brenn- 
ebene ist die Konvergenz der Integrale an der Grenze — oo erforderlich. 

Wenn wir nun voraussetzen, daB beide Brennebenen existieren, so k6nnen 
wir in formaler Weise nach nebenstehender 
Figur die entsprechenden -«Brennweiten » / 
aus den achsenparallel einfallenden Strah- 
len definieren. 

Diese Brennweiten errechnen sich aus 
den Formeln 


Fig. 1 
i= a(—o0) - f= o(+00) (9) Ubliche Definition des «Brennpunktes» als 
i 0 


0’ (2p) 


: Bild F, des unendlichfernen Dingpunktes, 
H, Hauptpunkt», F, H, «Brennweite». 
worin o(z) und o(z) die von links bzw. 

rechts achsenparallel einfallenden Strahlen bedeuten. Die Kenntnis von Brenn- 
punkten und Brennweiten gentigt aber, um zum Objekt G ein «Bild» B* nach 
der tiblichen Bildkonstruktion, wie in Fig. 2 angedeutet, zu bestimmen. 


Man erhalt dann eine « VergroBerung» B* nach der Linsengleichung 


g* 
6g ae A eg el mE (10) 
G £9 — Zr, ht 


1) Siehe z. B. W. GLASER, Z. Phys. 81, 637 (1933); 83, 104 (1933); M. Corre, Ann. Phys. [11], 
10, 333 (1938); und J. Laprume, Cah. Phys. 29 (1946). 
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Die auf Grund dieses Verfahrens bestimmte Abbildung wollen wir die zur 
elektronenoptischen Abbildung «zugeordnete Newtonsche Abbildung» nennen. 
Das so gewonnene «Bild» B* wird aber sowohl nach Lage und GréBe vom tat- 
siichlichen Bild, das nach der strengen Abbildungsgleichung (2) berechnet wird, 
verschieden sein. Berechnet man daher die Brennweite aus der experimentell 
bestimmten VergréBerung unter Benutzung der Gleichung 


44 — 2p, 


eer pa 


(11) 


Fig. 2 : 
Zur Festlegung der «zugeordneten» Newtonschen Abbildung. B ist das Bild auf Grund der tat- 
sdchlichen elektronenoptischen Abbildung, B* ist das mittels «Brennweiten» und «Hauptpunkten» 
konstruierte «zugeordnete» Bild. 


so wird man einen Wert erhalten, der im allgemeinen nicht mit dem aus (9) 
bestimmten tbereinstimmen wird. Nur in dem besonderen Fall, daB zugeord- 
netes Bild B* und tatsachliches Bild B nach Lage und GréBe zusammenfallen, 
ist es berechtigt, aus der gemessenen VergréBerung f (= f*) die Brennweite 
nach (11) zu berechnen. Stimmt zugeordnetes und tatséachliches Bild nicht 
liberein, hat die nach Fig. 1 bestimmte Brennweite keine physikalische Bedeu- 
tung und wird insbesonders nicht unmittelbar die VergréBerung bestimmen. 

Man erkennt aus Fig. 2, da elektronenoptische und zugeordnete New- 
tonsche Abbildung sicher immer dann zusammenfallen, wenn es sich um ein 
schwaches Abbildungsfeld handelt, d.h. wenn der Raum links von Fy, und 
rechts von /’, praktisch feldfrei ist. 

Die Aufgabe, alle starken Abbildungsfelder zu bestimmen, fiir welche die 
experimentelle Ermittlung der Brennweiten nach Gleichung (11) berechtigt ist 
und mit der nach Fig. 1 aus einem achsenparallel einfallenden Strahl gewon- 
nenen tbereinstimmt, wurde erstmalig von GLASER?) gestellt und dann in 
zwei weiteren Ver6ffentlichungen von GLASER und LAMMEL?) gelést. 


1) W. Graser, Z. Phys. 117, 285 (1941). 
2) W. GLAser und E, Lammet, Ann. Phys. 40, 367 (1941); Mh. Math. Phys. 40, 289 (1948). 
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_ In einem spateren Abschnitt dieser Arbeit sollen die magnetischen Abbil- 
dungsfelder dieser Klasse naher diskutiert werden. 


2. Oskulierende Newtonsche Abbildung 


Um die: Bedeutung der nach Fig. 1 ermittelten Brennweiten fiir die elek- 
tronenoptische Abbildung zu erkennen, wollen wir die allgemeine elektronen- 
optische Abbildungsgleichung (2), (3) in der Nahe zweier konjugierter Punkte 
untersuchen. Seien z) und 2, zwei nach (2) einander zugeordnete Punkte. Wir 
betrachten den Bildort z, + Az,, der zu einem benachbarten Dingort 2) + Az 
-gehort. 

Wir wollen unsere Betrachtungen zundchst auf rein magnetische Abbil- 
dungsfelder beschranken. Bei der Abbildung durch elektrische Linsen hat man 
es meist mit einem feldfreien Ding- und Bildraum zu tun, so da8 die Gesetze 
der zugeordneten Newtonschen Abbildung verwirklicht sind. Denn zur Vermei- 
dung von Stérungen des elektrischen Abbildungsfeldes mu8 man das Objekt 
auBerhalb desselben annehmen. Wie sich tibrigens zeigen wird, lassen sich fur 
starke elektrische Abbildungsfelder (Objekt im Feldbereich) die nachfolgenden 
Beziehungen auch nachweisen, wenn eine bestimmte zusaitzliche Bedingung 
erfiillt ist. 

Bezeichnen wir die Partikularlésungen von (5) mit w und w, so haben wir 


a= Py: o= Ow, 


Im rein magnetischen Abbildungsfeld konnen wir daher g mit w und o mit w 
identifizieren. 

| Bezeichnen wir die VergréBerung fiir die Objektlage 2 + Az,amit B, so”er- 
halten wir 


ai U(Zy + Azo) wz) +Azp) ’ 


u(% + Az) w(z+A%) p= wey eh) (12) 
UWS 


(Zo) w(zo) * 
Der Zusammenhang zwischen Dingweite A 2 und Bildweite 4z, (Abbildungs- 
- gleichung) ergibt sich zu 


u(z, + Az) w(% + Az) = ¥(% + Azy) w(z, + A2,) . (13) 
Entwickeln wir nach 42, und 4%, bis zu Gliedern 3. Ordnuag, so erhalten wir 


ff 
Uy Wy — Wy Uy + (uy Wy — Wy uo) Azo + (Uz Wy — Uo wy) Az, 
. 


il 
+ (wi wh — uf, w4) Azy Az, + = (uz Wo — Mo #1) (Aa) 


1 y 
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Die zweiten Ableitungen «” und w” kénnen wir auf Grund der Differential- 
gleichung (5), der « und w geniigen, eliminieren. Wegen der Abbildungsgleichung 


Uy Wy — Wy Uy = 0 


verschwindet daher das erste und letzte Glied. 
Wir schreiben zur Abkiirzung 


, 7 
Wy Up 2 ea ae | Wy Ug 
a= | , Pal ae, b= | , , | , a ) we ue (14) 
Wy Uy Wo Uo | 1 its 
und 
Wy Up W, Uy, |, : 
C = s EM a , ele (15) 
Wo Uo Wy 14; | } 


Die Konstanz des Ausdrucks (15) folgt in bekannter Weise, wenn wir die Diffe- — 
rentialgleichung (5) fiir « mit w und die fiir w mit ~ multiplizieren und dann 
subtrahieren. Zwischen den Konstanten (14) und (15) bestehen gewisse Zu- — 
sammenhange. Berechnen wir fC, indem wir das erste Glied von C mit 
B = w,/W) und das zweite Glied mit 6 = u,/u, multiplizieren, so erhalten wir 


[COE ae (16) 
Analog zeigt man 
pase Gu. (17) 
Daher ist 
b 
pra. (18) 
Die Abbildungsgleichung (13) lautet nunmehr 
b 
— Az 
Az, = - Ge8 set te (19) 
1+ — Az, 1+ AZ 
und die Vergr6Berung (12) 
1+ 4. Age 
ee er (20) 
1+ — Az, 
Ug 


Wenn wir aus der Abbildungsgleichung fiir A z, einsetzen und konsequent Glie- 
der héherer Ordnung vernachlassigen, bekommt Gleichung (20) die Gestalt 


B=B{1+ydAz}, (21) 
wobei wir zur Abkiirzung 


2 Uy UG ’ 
ea ea (22) 


Ug 


gesetzt haben. 
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Nach (15) gilt 


1 Wy Uy — C ; Uy W,— C 
Wo al Wy A ? 
| Uo Uy 
also ist 
* pes , 
WwW w U 
; cen Suse oe = 
Uo Cal Uo Uy 


Auf Grund der Abbildungsgleichung verschwindet der erste Ausdruck, und wir 
erhalten 


_— CC (ue uy Mm) _ 7 
Joel Avi female Cy 
so daB nach (17) 
c=-—ay (24) 
ist. Die Vergr6Berung ist daher 
B= 6{1— 2 Az} -—f—. (25) 
al =- Gy 4% 


Der Vergleich mit Abbildungsgleichung (19) zeigt, daB die «Tiefenscharfe» all- 
gemein in der untersuchten Naherung durch 


Az 
ae = BB (26) 


gegeben ist. Die gewonnenen Formeln (19) und (25) kénnen wir mit den ana- 
logen Formeln der gewohnlichen Lichtoptik vergleichen. Zwei konjugierte 
Punkte sind durch die Gleichungen 


21 — 2p, fe : 4 — Zp + AY i 
= = B= pe 0 
B h 29 — 2p, s h Zo — 2p, + Azo 2) 
verkniipft. Aus der zweiten Gleichung folgt 
re A B? Azo 
An = are ae : (28) 
1+ — Az, 
0 
und fiir die VergréBerung ergibt sich 
eee , : (29) 
1+ aes Az 
lo 


Wir erkennen, da8 unsere Formeln (19) und (25) die gleiche Gestalt haben wie 
die entsprechenden Formeln der Lichtoptik. Es zeigt sich weiterhin, daB wir, 
ohne zu einem Widerspruch zu kommen, die entsprechenden Ausdriicke gleich- 


ZAMP 1/24 
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setzen kénnen. Durch Gleichsetzen des Zahlers und des Nenners erhalten wir 


hh Se fo (30) 
und 
#. =<. (31) 


Durch diese beiden Gleichungen ist eine Newtonsche Abbildung definiert, welche 
wir die «oskulierende Newtonsche Abbildung» unserer elektronenoptischen’ 
Abbildung in den beiden konjugierten Punkten z, und z, nennen wollen. 

Fiir die Brennweiten der oskulierenden ew ECBO Abbildung ergibt sich 
nach (31), (17) und (14), (15) 


1 Wo uy — = wy U6 (32) 
fo Wo — WG Uo 

und nach (30) 
by Wy Wy — Wo Uy (33) 
ih W, Uy — Wy Uy” = : 


so daB 1//, aus 1/f) durch Vertauschung der Indizes 0 und 1 hervorgeht, wenn 
wir bedenken, daB die Wronski-Determinante (15) im Nenner bei dieser Ver- 
tauschung ungeandert bleibt. Da die Brechkraft durch den Quotienten zweier 
Determinanten aus zwei Partikularlésungen gegeben ist, bleibt sie ungeandert, 
wenn man zu einem neuen Paar von Partikularlésungen tibergeht. Das mu8 
nattirlich der Fall sein. 

Diese Brennweiten haben im Rahmen der allgemeinen elektronenoptischen 
Abbildung folgende Bedeutung: 

Ist v eine Elektronenbahn, die im zweiten Bezugspunkt z, achsenparallel ist, 
so kann sie mittels einer Konstanten durch 


ont rveenaede 


y = const (wu, — wu) (34) 
gegeben werden. Die Brennweite /, ist dann nach (32) durch 
n. r(2,) 
fy = et (35) 


gegeben. Wir betrachten ein Elektronenbiindel, das die zweite Bezugsebene x, 
senkrecht durchsetzt (Fig. 3). Man sieht sogleich, daB die Tangenten an die 
Elektronenbahnen dieses Biindels in der ersten Bezugsebene zy die Achse in 
ein und demselben Punkt schneiden. Diesen Punkt nennen wir den oskulieren- 
den Brennpunkt fy) der Abbildung. Die Entfernung der Ebene %» von diesem 
Brennpunkt bezeichnen wir mit Z,. Sie ist gegeben durch 


Foe oot tae (36) 
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Wenn wir fiir 7 aus (34) einsetzen, erkennen wir, daB die Konstante, welche den 
Achsenabstand der einzelnen Strahlen in der Ebene z, festlegt, herausfallt und 
Z, die Gestalt 
: Z, Do Baedie IC je 

Wo Uy — Wi Uy c 


(37) 


annimmt. Ganz analog ergibt sich fiir das achsenparallele Biindel in der Ebene 
Zt, als zweite «Brennpunktlage» 


W, Uy — Woy 
LZ, re Dy a Rants ame 2 
Wi Ug — Wo U4 c 


(38) 


Weiter gilt, daB die Schnittpunkte der obigen Bahntangenten in z mit den zu- 
_gehérigen achsenparallelen Bahntangenten in m7 eine achsensenkrechte Ebene 

bilden, welche wir die oskulierende Hauptebene nennen. Fiir ihre Entfernung 
-vom Brennpunkt finden wir nach Fig. 3 


ple (2) WyUy— Wu  C 
H y'(Z) WO UL — WL UG Ca 


(39) 


Fig. 3 
Zur Erlauterung der oskulierenden Newtonschen Abbildung fiir die beiden konjugierten Bezugs- 
ebenen 7% und 7}. 


Der Beweis dafiir, daB diese Tangentenschnittpunkte eine achsensenkrechte 
Ebene bilden, liegt darin, daB sich fiir alle Bahnen der gleiche Wert Zz, 
ergibt, und man erhalt nach (39) und (35) 


Ze. = fo. (40) 


Es schiene vielleicht natiirlicher, als oskulierenden Brennpunkt den Achsen- 
schnittpunkt S der Strahlen des Biindels zu definieren. Dieser Punkt hat jedoch 
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im allgemeinen Fall keine unmittelbare Beziehung zur Newtonschen Abbildung. 
Im Falle, daB der Punkt S bereits im feldfreien Raum liegt, fallt der oskulie- 
rende Brennpunkt im obigen Sinne mit S zusammen. Ganz analog ergibt sich 
die zweite Brennweite und Hauptebene. 

Unsere Abbildungsgleichung (19) kénnen wir in der Gestalt 


(Az + 24) (4% + Zo) = foh (41) 
schreiben. Durch Vergleich von (31) mit (37) ergibt sich 
Eom 2, (42) 
und fiir die VergroBerung erhalten wir nach (25) 
fe fo 
Bez 5 Ra ae (43) 


Wir erkennen daraus, daB die VergréBerung fiir 42) = /) — Zp eins wiirde. 
Az = fy — Zp gibt uns aber gerade die Entfernung der oskulierenden Haupt- 
ebene von der ersten Bezugsebene zp. Damit haben wir fiir unsere oskulierenden 
Haupt- und Brennpunkte die gleichen Eigenschaften nachgewiesen wie in der 
Lichtoptik, und es kann daher auch die bekannte Listingsche Bildkonstruktion 
angewendet werden. 

Wir wollen nochmals darauf hinweisen, daB strenggenommen jedes konju- 
gierte Punktepaar seine eigenen oskulierenden Kardinalelemente besitzt!). Nur 
in den Fallen der schwachen Abbildungsfelder oder der Felder mit Newtonscher 
Abbildungsgleichung fallen diese oskulierenden Brenn- und Hauptpunkte fiir 
alle konjugierten Punktepaare zusammen. 

Die experimentelle Bestimmung dieser oskulierenden Brennweite kann 
strenggenommen — d.h. wenn das Feld nicht als Feld mit Newtonscher Ab- 
bildungsgleichung aufgefaBt werden kann — in exakter Weise nur nach der aus 
(43) folgenden Gleichung 


1 1S yal 1 

Teale “ 
durchgeftthrt werden. Die Objektverschiebung kénnte in der Art, wie dies 
E. Ruska?) gemacht hat, durch besondere MeBsiitze verschieden hoher Blen- 
denfassungen erreicht werden. Die neue Scharfstellung des Bildes hatte aller- 


dings nicht durch eine Abgleichung des Objektivstromes zu erfolgen, sondern 
durch eine entsprechende Verschiebung des Bildschirmes. Infolge von (26) 


1) Die oskulierenden Kardinalelemente fiir nichtnewtonsche Abbildungsfelder, zum Beispiel fiir 
die vom erweiterten Glockenfeldtypus B, = By/[1 + (z/a)2]# (u> 1), sollen in einer folgenden, 
gemeinsam mit F. Lenz verfaBten Arbeit explizit bestimmt werden. 


a) E. Ruska, Uber den Bau und die Bemessung von Polschuhlinsen fiir hochauflésende Elektro- 
nenmikroskope, Arch. Elektrotechn. 38, 102 (1944) 
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-wiirde bei einer VergréBerung 6 ~ 100 eine 10000fache Bildschirmverschie- 
bung notwendig sein, was bei der geringsten von Ruska erreichten Objektver- 
schiebung von 0,1 mm bereits eine Bildschirmverschiebung von 1 m erfordern 
wiirde. Infolge der hohen Tiefenscharfe wird natiirlich eine geringere Verschie- 
bung des Bildes geniigen. Vielleicht ist es bei einer kleineren Objektverschie- 
bung, die natiirlich anzustreben ware, sogar méglich, den Auffangschirm an 
seiner Stelle zu belassen. 
Fir den bei den Anwendungen besonders wichtigen Fall hoher VergréBerung 
k6nnen wir als die beiden konjugierten Bezugsebenen zwei Ebenen betrachten, 
die mit der im Sinne von Abschnitt 1 definierten Brennebene und der zugeord- 
neten unendlich fernen Ebene zusammenfallen, wobei wir aber noch immer 
mit einem endlichen Wert von f, zum Beispiel 100, rechnen diirfen. Der zuge- 
horige oskulierende Brennpunkt fallt in diesem Fall mit dem Brennpunkt im 
Sinne von Abschnitt 1 zusammen. Auch die Brennweite stimmt dann mit der 
friiher gebrachten (Fig. 1) tiberein. 
Wenn wir nun auch die Abbildung in einem starken elektrischen Abbil- 
_dungsfeld in Betracht ziehen, dann kénnen wir genau wie oben zeigen, daB die 
Zuordnung zwischen Ding- und Bildpunkt in der Nahe zweier konjugierter 
Punkte in derselben Approximation wie oben durch die Gleichung 


bi re (45) 


gegeben ist, wobei a,b und c mit den partikulaéren Integralen w und w der 
Differentialgleichung (5) nach (14) erklart sind. Wenn wir rein formal eine 


«VergréBerung» B und B durch 
+ Az;) u(z, + Az») 


NR rene. pain aay 
tae eer wey Tue a) 


einfiihren, haben wir dieselben Verhiltnisse wie friiher, und B ist daher durch 
jsp eos ee (47) 


gegeben. Die tatsdchlichen VergréBerungen B und f nach (3) stehen mit B und 
B in folgendem Zusammenhang 


7 Both, tae ey 
Ba ol t+4a) ees B- ; Ae 
0(Z + 425) OD, nore 


4 @, 
A lee 
o(%) _ D 7 | 
B ~~ (29) “ V2 B- 


(48) 
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Durch Einsetzen von (45) und (47) ergibt sich 


B a 
APA 1 & at a 
le+oe 4 OQ, a : 


Wenn wir in die Formeln (16), (17) und (18), welche fiir 8 gelten, die Vergr6Be- 
rung einfiihren, erhalten wir 


t 
£ = = ee 
O, - _ >}: Die pe 2/ = 2 50) 
pS cas; p/fianmc: pg =2 (50) 
Damit kénnen wir unsere Abbildungsgleichung (45) in der Gestalt schreiben _ 


B? oe Az, ] 
Az, == u . (51) q 
teeteee Az : 

a 4 

Vergleichen wir diese Gleichungen (49) und (51) mit den entsprechenden New- : 


tonschen (28), (29), dann sehen wir, da unsere Gleichungen eine oskulierende — 
Newtonsche Abbildung darstellen, wenn folgende Bedingungen erfiillt sind: 


@ 

al oa 
ee ee oe 

ne : (ae ao Le, 4: Ge la 


(53). 


Die erste Gleichung stellt die bekannte Beziehung zwischen den Brennweiten 
eines optischen Systems der gewohnlichen Lichtoptik dar. In der Elektronen- 
optik wird sie gewéhnlich nur unter der Voraussetzung von feldfreiem Ding- 
und Bildraum hergeleitet, wo sie unmittelbar aus der Konstanz von 


V'® (9 0’ — 90) (54) 


gefolgert werden kann. Man sieht, daB sie fiir die beiden oskulierenden Brenn- 
weiten allgemein gilt. Allerdings zeigt die zweite Bedingung (53), daB nicht 
jedes elektrische Abbildungsfeld in allen konjugierten Punktepaaren eine osku- 
lierende Newtonsche Abbildung besitzt. Es muB zwischen den Feldern in zwei 
zugeordneten Punkten die Bedingung 

D5 Do bd 

oor ee 0 (55) 


erfiillt sein. Mit (50) folgt hieraus 


D5 9 Dy 
ie ANE (56) 
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Nach (46) und (48) ist 


; Re Os =\% eee 
o P, Wo P, Uy” 
Mit Hilfe der Differentialgleichung (5) [fiir B,(z) = 0] koénnen wir daher die 
Bedingung (56) fiir die Existenz einer oskulierenden Newtonschen Abbildung 
im starken elektrischen Abbildungsfeld durch 


(w,\4- (m\t ) 5 
($)'= (y= (57) 
ausdriicken'). Eine genauere Diskussion dieser Bedingung und die explizite 
Bestimmung derartiger Felder wollen wir einer spateren Arbeit vorbehalten. 

Um die geometrische Bedeutung unserer optischen Konstanten im starken 
elektrischen Abbildungsfeld zu erkennen, miissen wir die fiir die Abbildung 
kennzeichnenden Determinanten a, 0, c und C durch die wirklichen Elektronen- 
bahnen 9 und o ausdriicken. Wir definieren dazu in analoger Weise 


oy ede Piel ky Bale ot: (58) 
A @1 1% @o| oe e@| 


Es zeigt sich, daB gerade auf Grund der Bedingung (55) die Beziehungen gelten 


CO ee a ee (59) 
0D; 


-z __|% Qo| rm _ | @1| 
oi mn Ca Ab ey 
so ist allgemein 
ra Cc a Cc 
= —— und C, = —= 61 
0 ) ®, 8 Vo, ( ) 
Fiir die reziproken Brennweiten findet man daher 
1 /@ 1 |/@ 
1 c 0 Cc 
Le es a2) a / eae 62 
lle xa js ®, C oS 
und wenn wir c, C in 1/f, durch c, C, ersetzen 
1 & _ Hat~ 0% _ 1 (40) 63) 
to Cs 0; OY — 191 v(%) ” 


womit gezeigt ist, daB die Brennweiten dieselbe geometrische Bedeutung haben 
wie im rein magnetischen Fall. 


1) Die Bedingung (55) fur die Existenz einer oskulierenden Newtonschen Abbildung ist trivialer- 
weise erfiillt, wenn ®,/D, = 0 und O;/M, = 0 ist, dh. Ding und Bild in feldfreien Gebieten liegen. 
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Die in gleicher Weise wie friiher definierten Brennpunktkoordinaten Z) un 
~ Z, sind durch 


(64) 


a1] oO 


z= ud Zo 


gegeben. Auf Grund der Formeln (59) erhalten wir auch hier wieder die Glei- 
chungen (37), (38), woraus die elektronenoptische Bedeutung der eingefihrten 
Determinantenverhiltnisse a/c und b/c hervorgeht. 

Wir sehen also, daB die Bedingung (55), (56) oder (57) das Bestehen aller 
Beziehungen der oskulierenden Newtonschen Abbildung sicherstellt. Gegen-— 
tiber den Verhdltnissen im rein magnetischen Abbildungsfeld, wo eine oskulie-_ 
rende Newtonsche Abbildung fiir zwei beliebige konjugierte Punkte immer 
vcrhanden ist, unterscheidet sich das elektrische Feld darin, daB dazu noch 
die besondere Bedingung (57) erfiillt sein muB. : 

Verlangen wir, daB die Abbildung im rein magnetischen Abbildungsfeld 
durch eine oskulierende Newtonsche Abbildung in héherer Ordnung approxi- | 
miert wird, dann fiihrt das, wie im Anhang II gezeigt wird, zu der Bedingung | 


(=)'- (2 y Qo Biz) (65) | 


Wo Uo - Q, i Biz) ’ 


welche bemerkenswerterweise mit der Bedingung (57) des rein elektrischen 
Feldes fiir eine um eine Stufe niederere Approximation tibereinstimmt. 


Anhang I 


Bedingungen fiir die Existenz von geradlinigen Asymptoten’) 


VoraussetzungsgemaB soll fiir den unendlich fernen Punkt z +o die Tan- 
gente 


@ — e(2) = @'(2) (z— 2) (I, 1) 
existieren, also folgt 
lim o’(z) = a; lim [@(z) -- o’(z) z] = 5, (TZ) 


wobei a und 6 bestimmte Konstanten sind. Ausgedriickt in R lauten diese 
Bedingungen nach (4) 


lim "14 (Rg) — 2% RG) = a (I, 3) 
lim @-18 [Re) (1 + 7-5-2) — B®) z| = 8. (I, 4) 


1) Den Beweis verdanken wir Herrn P. Scuiskr, 
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Da aus physikalischen Griinden ® weder Unendlich noch Null werden kann, 
wird R zugleich mit @ nach (4) unendlich. Da R’ stets endlich bleibt, mu8 also 
gleichzeitig mit (I, 2) nach (I, 3) 


lim R’'(z}= A und lim®R=A (i>) 
sein. Aus der zweiten Gleichung schlieBt man, daB ®’ im Unendlichen ver- 
_schwindet. Wir wollen nun voraussetzen, daB es starker als 1/z verschwindet. 


Dann lautet (I, 4), da ® gegen eine Konstante geht, 


lim (R— R’2)=B. (I, 6) 


200 


Diese Bedingung ist aber zusammen mit (I, 5) dafiir kennzeichnend, dab die 
Lésungen der Differentialgleichung (5) geradlinige Asymptoten besitzen, wie 
man durch Vergleich mit (I, 2) erkennt. 

Wir multiplizieren (5) mit z und erhalten durch partielle Integration 


R-Riz=[QRedz. (I, 7) 


Bilden wir den Limes (I, 6), so erhalten wir fiir ein geniigend groBes « 
fee R’ 2) —-B = lim [QRad: 


z—>0O a—>0o 


e 
p 


=| QRzdz + | Q(B- R'z)edz=J+ B[Qzdz- A [Qztdz, 


wobei J eine Konstante bedeutet. 
Fiir die Existenz einer achsenparallelen Asymptote geniigt also die Kon- 
vergenz des Integrals 


Jord. (I, 8) 


und fiir die Existenz beliebig geneigter Asymptoten ist die Konvergenz des 
Integrals 


co 


/Q dz (I, 9) 
erforderlich. . 
In analoger Weise, wie in dem Buch Gruppentheorie und Quantenmechantk, 
S. 63, von H. Wevt geschlossen wird, kann man zeigen, daB diese Bedingungen 
auch hinreichend sind, wobei man nur in unserem Falle als Ausgangslésung des 
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in dem erwahnten Buche angefiihrten sukzessiven Approximationsverfahrens 


R,=Az+B (I, 10) 
zu setzen hat. 


Anhang II 


Oskulierende Newtonsche Abbildung im rein magnetischen 
Abbildungsfeld in héherer Naherung 


Um zu den Formeln der oskulierenden Newtonschen Abbildung in héherer 
Naherung zu gelangen, betrachten wir die Bildkoordinate und die zugeh6rige 
VergréBerung als Funktionen der Dingkoordinate und entwickeln beide in der 
Umgebung des Bezugspunktes z, in eine Taylor-Reihe. Aus der allgemeinen Ab- 
bildungsgleichung im rein magnetischen Abbildungsfeld 

B (zo) = w(lZ (Zo) } ae ul (Zo) | (II, 1) 


w (Zo) ~ u(Z) 


worin 2, als Funktion von z, aufgefaBt wird, erhalten wir durch Differentiation 
nach Z) auf Grund von (14) und (18) 


dz, 
rE ege un 
und 
dp C pe 
Fr epee on le es 


C und ¢ sind die in (14) und (15) eingefiihrten Determinanten. Aus (II, 2) und 
(II, 3) folgt unmittelbar 
Oe deed pu (II, 4) 
dz (Gags ; 


Bildet man die dritte Ableitung von z,, so ergibt sich 
d8z b 
Te ~ ~8 BSE E+ 2h Oe - 28S S2, LS) 
also 
Giz 


is a pat — pe Bt), (IT, 6) 
In den beiden letzten Formeln wurden die Differentialgleichung (5) und die 


Formeln (16), (17) und (18) benutzt. 
Ganz analog erhalt man 


ger = 28°(G) + BQo(1— Bt &). (I, 7) 


7 - 
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Die Abbildungsgleichung und die Vergré8erung sind daher in dieser Naherung 
zegeben durch 


A 2, = %(% + A%) — %(20) = B? 42 {1 amie ra Az 
, i o\2 1 Q i (II, 8) 
+ [o(Zh+ 4 proalt— pe SV] rote} | 


B-p{i—£64x + [a(E)+ + Oo(t— Bt $)| Ga)?+-}. 9) 


Durch Vergleich mit den entsprechenden Reihenentwicklungen der Newton- 
schen Abbildungsgleichungen (28) und (29) erkennt man, daB die elektronen- 
optische Abbildung im Magnetfeld eine oskulierende Newtonsche Abbildung 
bis inklusive Glieder 3. Ordnung darstellt, wenn die Bedingung 


Bs =(Sy- (4)- Qo _ B:{20) (IT, 10) 


Wg / Mg} 0, BF (41) 


erfiillt ist oder wenn Q, und Q, verschwinden. 


Résumé 


La notion de distance focale est étudiée pour les lentilles électroniques. La 

grandeur de la diminution de champ, qui est nécessaire pour l’existence de foyers 
et de distances focales, est déterminée. La définition habituelle de la distance 
focale ne caractérise pas en général la formation de l’image dans le cas de lentilles 
électroniques fortes. La relation entre la position de l’objet et celle de Vimage, de 
méme que le grandissement linéaire sont rendus, en optique électronique, par ses 
fonctions plus compliquées qu’en optique de la lumiére. Ces fonctions sont appro- 
chées, dans le voisinage de deux points conjugués par l’«équation newtonienne | 
osculante pour la fonction de l'image» jusqu’aux termes du quatriéme ordre 
inclusivement. A chaque paire de points conjugués correspondent ainsi les foyers, 
points principaux et distances focales résultant de l’«équation newtonienne oscu- 
lante». Si les éléments cardinaux osculateurs sont indépendants de la paire de 
points conjugués choisie, ils caractérisent & eux seuls la formation de Vimage et 
sont identiques aux grandeurs définies de la facon habituelle. Nous avons appelé 
de tels champs, pour lesquels 1’équation de formation de l'image de l’optique 
ordinaire est strictement valable, champs de représentation newtoniens. 

La seconde partie de ce travail est consacrée a l’étude de ces champs forts. 
Des exemples de tels champs, qui approchent des champs existant réellement 
dans dés lentilles électroniques, sont donnés. Si l’on veut conserver a la notion de 
distance focale un sens physique pour n’importe quels grandissements, il faut 
approcher aussi bien que possible le champ empirique par un champ de la classe 
newtonienne. On indique un procédé pour le faire. Enfin on examine des méthodes 
expérimentales, qui permettent de déterminer les distances et les points focaux 
de ces champs approchés ayant une formation d’image newtonienne. 


(Eingegangen: 28. 2. 1950.) 
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Fligelschwingungsformen 
in ebener kompressibler Potentialstromung 


Von Nikoraus Rott, Ziirich*) 


: 
Wenn bei einem schwingungsfahigen System*von  Freiheitsgraden nur 
elastische Krafte und Tragheitskrafte auftreten, so ergeben sich Trecmeaeell 
der n Hauptschwingungsformen als Lésung der sogenannten Sakulargleichung, ~ 
einer algebraischen Gleichung m-ten Grades mit reellen Koeffizienten, die zudem 
die bekannte Besonderheit hat, 1 reelle Lésungen zu besitzen. Beim Flatter-_ 
problem wird die Aufgabe verwickelter durch das Auftreten von Luftkraften, — 
wodurch die Bestimmungsgleichung fiir die Frequenz komplexe Koeffizienten — 
erhalt. Es ergeben sich im allgemeinen auch komplexe Frequenzen als Lesude 
so daB je nach dem Vorzeichen des Imaginarteiles der Frequenz Dampfung 
oder Anfachung eintreten kann. Die komplexen Koeffizienten hangen dabei 
noch ab von der Fluggeschwindigkeit. Bei einer gewissen sogenannten «kri- 
tischen» Geschwindigkeit kann nun der Fall eintreten, daB eine Frequenz rein 
reell wird, die entsprechende Schwingungsform also weder gedémpft noch an- 
‘gefacht ist. Die Bestimmung der kritischen Geschwindigkeit des Fliigelflatterns - 
bedeutet algebraisch die Lésung der folgenden Aufgabe: eine Gleichung ist 
gegeben mit komplexen Koeffizienten, die von einem Parameter abhangen; 
gesucht ist derjenige Parameterwert, fiir den eine rein reelle Lésung existiert. 

Bei der Durchfiihrung dieser Aufgabe zeigt es sich, daB die Rechnungen 
schon bei zwei Freiheitsgraden ziemlich umfangreich werden und daB insbe- 
sondere der Zusammenhang zwischen den Ausgangsgr6Ben und der kritischen 
Geschwindigkeit uniibersichtlich wird. 

In den folgenden Ausfiihrungen wird untersucht, welche Schwingungsformen 
bei der kritischen Geschwindigkeit auftreten kénnen. Es zeigt sich, daB die 
moglichen Schwingungsformen durch die Art der Luftkrafte allein eingeschrankt 
sind. Die méglichen Schwiugungsformen im kritischen Falle sind diejenigen, 
bei denen die Luftkrafte weder eine Arbeit leisten noch Arbeit aufnehmen, das 
heiBt die weder angefacht noch gedémpft sind. Die elastischen und Tragheits- 
krafte leisten ja infolge ihres konservativen Charakters insgesamt bei periodi- 
scher Bewegung keine Arbeit. Die wirkliche Schwingungsform ergibt sich unter 
den méglichen natiirlich erst bei Lésung der Bewegungsgleichung, das heiBt 
der eingangs beschriebenen Aufgabe. Aber selbst durch Einschrankungen kann 
man den Umfang dieser Aufgabe vermindern, was besonders willkommen sein 


1) Institut fir Aerodynamik, ETH. 
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diirfte bei der Erhéhung des Rechenaufwandes, den die Beriicksichtigung der 
Kompressibilitat mit sich bringt. 
Die Methode, mit Hilfe der Arbeitsberechnung die méglichen Flatterformen 
‘zu untersuchen, ist schon von GREIDANUS, DuNCAN (inkompressible Strémung) 
und Barton (Uberschallstrémung) angewendet worden. Anscheinend hat je- 
doch diese bemerkenswerte Betrachtungsweise noch keine sehr groBe Verbrei- 
tung gefunden, so daB sie im folgenden ausfiihrlich dargestellt werden soll. 
Die Ergebnisse werden diskutiert fiir kleine Schwingungen in ebener Potential- 
str6mung im ganzen Mach-Bereich von Obisoo, wobei die Luftkraftkoeffizienten 
der vorhandenen Literatur entnommen wurden, mit Ausnahme des Sonderfalles 
_M = 1, fiir den auf Grund eigener Rechnungen Tabellen der Beiwerte angefiihrt 
sind. — Die Diskussion fiihrt zu einer von der Mach-Zahl abhangigen oberen 
Grenze der reduzierten Frequenz, tiberhalb welcher keine angefachten Schwin- 
gungen mit zwei Freiheitsgraden méglich sind, und zu einer niedrigeren, im 
Uberschallgebiet zum Teil schon bekannten Grenze fiir Schwingungen mit nur 
einem Freiheitsgrad. Die Betrachtungsweise vermittelt auch ein anschauliches 
Bild des Flattervorganges und der MaBnahmen, die zu seiner Verhiitung ge- 
troffen werden kénnen. 


1. Die Luftkraftkoeffizienten und ihre Transformation 


Im folgenden beschranken wir uns auf das Flatterproblem mit zwei Frei- 
heitsgraden in ebener Stromung. Ein Profil der Tiefe / sei mit der Geschwindig- 
keit U angeblasen. Zur Kennzeichnung seiner momentanen Lage wird ein Be- 
zugspunkt auf der Profilachse gewahlt (P in Fig. 1); die Lagekoordinaten sind 
der Ausschlag y des Bezugspunktes und der Anstellwinkel « der Nullauftriebs- 
achse. 


Lagekoordinaten und Luftkrafte am. Profil. 


Unter dem Einflu8 von Luftkraften, Tragheitskraften und von elastischen 
Riickstellkraften, welche den Ausschlagen y und « entgegenwirken, kann nun 
der Fliigel Schwingungen mit zwei Freiheitsgraden ausfiihren: Schlag- und 
Drehschwingungen. Es werden dann y und « harmonische Funktionen der Zeit 
mit einer Kreisfrequenz v, komplex geschrieben: 


aye, a Oye" (1) 
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(Um Gleichdimensionalitat der Koordinaten zu haben, beniitzen wir spater 
yl.) Die Amplituden miissen im allgemeinen zur Kennzeichnung der Phasen- 
lage komplex angenommen werden. 

Die Luftkrafte kénnen nun fiir die durch (1) gegebene Bewegung des Profils 
bei kleinen Ausschlagen in einer Potentialstrémung berechnet werden, sowohl 


inkompressibel als auch fiir kompressible Medien. Als Ergebnis einer solchen — 


Rechnung kennt man den Auftrieb 2{ und das Moment St um den Bezugs- 
punkt. Beide werden ebenfalls harmonische Funktionen der Zeit mit dem 
gleichen Zeitfaktor e'”'; wir beniitzen die entsprechenden Beiwerte: 

W ivt mM tvt 


oe TTD ae ek oe eRe ae © (2) 


Vorlaufig brauchen wir als einziges Ergebnis, da sich fiir kleine Ausschlage 
zwischen ¥, und a% einerseits und c,, und c,,, anderseits ein linearer Zusammen- 
hang ergibt: 


Vo 


y , 
CA, an Ray - = Ra Ko 5 Cong - Bins 1 a Reine Xo. (3) 


Die k,, usw. sind komplexe Koeffizienten, da im allgemeinen zwischen yp, %% 
und ¢q,, Cm, Phasenverschiebungen bestehen; in Real- und Imaginarteil zerlegt, 
sollen die Koeffizienten in der Form geschrieben werden: 


Ray Pa Rare +4 Ravi (3a) 


usw. Man braucht also insgesamt acht reelle GroBen, «Luftkraftkoeffizienten », 
um den durch die Gleichungen (3) gegebenen Zusammenhang zu beschreiben. 

Fur ebene kompressible Potentialstr6mung ergeben sich die Luftkraftkoeffi- 
zienten in Abhangigkeit von zwei Parametern: 

1. die Machsche Zahl M der Anstrémgeschwindigkeit, 

2. die reduzierte Frequenz‘) 


yl 
(Oy 7 ; (4) 


w laBt sich in bekannter Weise auch als geometrischer Ahnlichkeitsparameter 
deuten. Betrachten wir das Profil in ruhender Luft mit der Geschwindigkeit U 
bewegt, so ist der wahrend einer vollen Schwingung zuriickgelegte Weg L gleich 


Cao : 
Vv 
damit wird auch 
l 
== 2 ia ; (4a) 


ausgedriickt durch das Verhialtnis der Profiltiefe 7 zur «Wellenlange» L. 


1) In einem Teil der Literatur wird die mit der halben Tiefe berechnete reduzierte Frequenz 
bentitzt: wo =v 1/2 U. 
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SchlieBlich hangen die Werte der Luftkraftkoeffizienten auch davon ab, 
welchen Bezugspunkt wir gewahlt haben. Wir befreien uns von dieser trivialen 
‘Abhingigkeit durch Angabe des Transformationsgesetzes beim Ubergang auf 
einen neuen Bezugspunkt P’ (Fig. 1), der ebenfalls auf der Profilachse liegen 
soll, jedoch gegeniiber P eine Riicklage hat vom Betrage 


PP’ =x. 


Alle auf den neuen Bezugspunkt bezogenen GréBen seien durch einen Strich 
gekennzeichnet. Ein Blick auf Fig. 1 bestatigt das folgende Transformations- 
gesetz fiir die Koordinaten (bei kleinem «) : 


y= y—ula, aw =o (5) 
und fiir die Krafte und Momente: 
PS, MS Mt oe 
oder mit Verwendung der Beiwerte: 
PO Cea GE He, (6) 


Durch Einsetzen von (5) und (6) in (3) erhalt man fiir die entsprechenden 
Amplituden 


C= Ray 2% + (Ray + % Fay) % 
Cn. = (Rin y ns Ray) * ve [Rana + %# (Raw - ny) + 2? kay] oa , 


wodurch in Analogie zu (3) neue Luftkraftkoeffizienten hk,’ usw. definiert sind, 
fiir die sich das folgende Transformationsschema ergibt: 


Ruy oa Ray , | 
Raw = Raw i % ky yp? 
. (7) 
tes rer ae 
Rena an Rime +H (Raa sa Remy) + x? Ray , 


Die Gleichungen gelten sowohl fiir die Real- als auch fiir die Imaginarteile. 


2. Das komplexe Amplitudenverhaltnis 


Die Beniitzung unndtig vieler komplexer GréBen wird vermieden, wenn 
man das Amplitudenverhaltnis 


ba (8) 
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in die Gleichungen (3) einfiihrt; ¢ selbst ist im allgemeinen komplex : 
CaS ee (8a) 
(, 7 reell). Schreiben wir nun (3) in der Form 
Co, = %0 Fay 6+ Fae) > Sy = %0 Amy SF Fame) (9) 


so bedeutet es keine Einschrankung der Allgemeinheit, wenn wir die Phasen- 
lage so normieren, daB a, rein reell wird. Dies soll fortan vorausgesetzt werden. 
Zur iibersichtlicheren Gestaltung der spateren Rechnung setzen wir in (9) 


Ray 6 of Raa a h, , Rny C mi Renee =a Vy , (10) 
wobei wieder 
hg = Nar + tha: him, = Vine + thmé (10a) 


geschrieben wird. Durch Identifikation von Real- und Imaginarteil in den 
Gleichungen (10) erhalt man mit den Bezeichnungen von (3a), (8a) und (10a): 


Nay = Rays om Ra yi H+ ee , 
hy a Fag g a ie Ui a Rawi , 
Mune = Biggns = my i n+ keer 
Nin; = Rm yi é ale Lee Ul a Rn wi c 


(11) 


Real- und Imaginarteil des komplexen Amplitudenverhaltnisses gestatten 
nun eine einfache geometrische Deutung, die schon von KASSNER und FINGADO 
[1]') angegeben wurde. Wir bemerken zunachst, daB nach (5) fiir € das folgende 
Transformationsgesetz gilt: 


C=C-—x (12) 

oder 
a (12a) 
y= y. (12b) 


Der Imaginarteil 7 ist also eine Invariante. Wenn wir nun einen neuen Bezugs- 
punkt wahlen derart, daB seine (auf die Lange / bezogene) Riicklage x = & ist, 
so wird wegen (12a) das Amplitudenverhaltnis rein imagindr. Diesen neuen 
Bezugspunkt nennen wir nach KAssNER und FINGADO den « Knotenpunkt». 
Fiir den Ausschlag yx des Knotenpunktes und den Winkel O&K = % (inva- 
riant) gilt also 
ie 
XK J 


= 17 


oder .. 
a) a vy tvt 
Ves toalyne =a 


) Die Zahlen in eckigen Klammern verweisen auf das Literaturverzeichnis auf S. 410. 


Saget gaa 


me hdenare 
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Im Knotenpunkt besteht zwischen Ausschlag und Verdrehung eine Phasen- 
verschiebung von 90°, reell geschrieben: 


Vr =—lHsinvi, x=) cosvi. 


Betrachten wir nun die absolute Bahn des Knotenpunktes vom Standpunkt 
der ruhenden Luft durch Einfiithrung des Weges s = Ut, den das Profil wah- 
rend der Teit ¢ zuriicklegt. Es ist mit Verwendung der reduzierten Frequenz 
nach (4) 


Ye = al nsino —. (13) 
Die Anstellung des Profils ist analog in Funktion von s gegeben durch 
= ag cos@ =. (14) 
Die Neigung # der Knotenpunktbahn ist nun gleich 
p= — = —Xy) W 4 COSW -. (15) 
Aus (14) und (15) folgt (mit 6, = Amplitude von £) 
= a =o 7) (16) 


das heiBt Neigung des Profils und Neigung der absoluten Bahn haben fiir den 
Knotenpunkt die gleiche Phase; das Verhaltnis der Amplituden dieser Nei- 
gungen ist durch 7 und w nach (16) gegeben. 

Durch die drei GroBen w, & und 7 ist eine Schwingungsform (bis auf einen 
Proportionalitatsfaktor x) eindeutig festgelegt. Fig. 2a zeigt folgendes Beispiel: 
Bezugspunkt sei die Profilmitte, und es ist o=)9§ =—0,5 und 7==0,5 
angenommen. Aufgezeichnet ist die absolute Bahn des Knotenpunktes, der 
wegen = —0,5 auf der Profilvorderkante liegt. Aus w folgt das Verhaltnis 
der Profiltiefe zur Wellenlange, fiir # =1: /= L/2a. Die Neigung des Profils 
ist tiberall proportional zur Knotenbahnneigung, sie ist also Null, wenn yx 
maximal ist, und hat den gréBten Wert fiir yx = 0 . Das Verhaltnis der Ampli- 
tuden, das im Punkt yx = 0 abgelesen werden kann, ist in unserem Beispiel 
nach (16) Bo 


Xp 


=-—wn=0,5. 


Ein anderes Beispiel sei w = 1, = —0,5, 7 =0. Wegen 7 = 0 ist nach 
(13) yx =0 und 6 =0. Man sieht, daB in diesem Sonderfall der Knotenpunkt 
zum Drehpunkt einer reinen Torsionsschwingung wird (Fig. 26). Reine Schlag- 
schwingung ist der «Sonderfall des Sonderfalles» mit dem Drehpunkt im Un- 
endlichen, also7 =0,€=+o. 


ZAMP 1/25 


386 NrkoLaus Rott ca 


Um die Bedeutung von positiven 7-Werten zu zeigen, ist noch ein Beispie 
mit a =1,&=—0,5,7=0,5 in Fig. 2e dargestellt. 


Z 


Fig. 2 
Fliigelschwingungsformen, dargestellt durch die Lage des Profils an verschiedenen Stellen der 
Knotenpunktbahn. Kennzeichnende GréBen (Bezugspunkt: Profilmitte) : 
a) o=1; €=—035; 7 =—0,5 
b)@=1; €=—05; 7= 0 
c)@=1; €=—05;7= 05 


3. Die Arbeitsleistung der Luftkrafte 


Die Frage, wann die Luftkrafte und Momente dampfend oder anfachend 
wirken, laBt sich nicht durch einen Blick auf die Gleichungen (3) beantworten. 
Die Antwort hangt ab von der Phasenlage von c,, und c,,, gegen Vo und %, und 
diese wird untibersichtlich, wenn die k,,, usw. komplex sind und die Ausschlage 
ein komplexes Verhaltnis haben. Es muB deshalb, wie schon einleitend bemerkt 
wurde, die Arbeitsleistung der Luftkrafte ausgerechnet werden. Die momentane 
Arbeitsleistung ist allgemein: 


v — (Kraft) - (Geschwindigkeit des Aufpunktes) 
+ (Moment um den Aufpunkt) - (Winkelgeschwindigkeit) 
Wy +- Ma 


oder mit den Luftkraftbeiwerten nach (2): 


1 eee ; 
2= > 0 UNF (c, = ea tl. (17) 
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Dabei ist nach (1) und (8) 
; | er ere a=tvae”. (18) 
Die c, und c,, seien in der durch (9) und (10) definierten Form beniitzt: 

Ce Ag hae”, Cy = agh,e”. ' (19) 


Bei der Ausrechnung von 2 mu8 man bekanntlich wieder auf die reelle 
Schreibweise iibergehen, das hei8t von den Ausdriicken (18) und (19) ist der 
Realteil zu nehmen unter Beriicksichtigung der Normierung, da8 a» rein reell 
ist. Mit den Bezeichnungen von (8a) und (10a) wird reell geschrieben: 

+ = —7 % (ny cosvt + E sin?) , a= —Vvasinve, (18a) 
Ca = Xo (Agr COSVE— gz Sin vt), Cy = % mr COSY E — hy; Sin yt) - (19a) 


Man erhalt & nach Einsetzen von (18a) und (19a) in (17). Sehr einfach wird der 


Ausdruck fiir die mittlere Leistung 2, der sich unter Beriicksichtigung der 
Mittelwerte 


5 il aon ee 
cos*#t= sin? vi =—, cosytsiny t = 0 


ergibt: 


~ i 
= = UP y 0G (— har + § hat + Ami) - 


Wir definieren nun einen «Leistungsbeiwert» A gemaB 


pi 


outrage — OM har + F htas + bmi) - (20) 


y = 
Durch Einsetzen von /,, usw. nach (11) in (20) ergibt sich nach leichter Um- 
formung: 


A = wo { Ray; ( r 7”) a Rnyi a Rai) g + Page oa fe, 3p) Ui + heat (21) 


Die Gleichung (21) gestattet die Berechnung des Leistungsbeiwertes fiir eine 
gegebene Schwingungsform, gekennzeichnet durch w, € und 7, wenn die Luft- 
kraftkoeffizienten bekannt sind, die selber von w und der Mach-Zahl abhangen. 
Explizit ersichtlich ist jedoch in (21) die Abhingigkeit des Leistungsbeiwertes 
von é und ny: es ist 4 eine quadratische Form in & und 7. 

Ist A positiv, so wird von den Luftkraften wahrend einer vollen Schwingung 
im Mittel Arbeit geleistet, die entsprechende Schwingungsform ist also ange- 
facht, wihrend umgekehrt negatives A Dampfung bedeutet. Im besonderen 
Falle A= 0 tritt weder Dampfung noch Anfachung. ein. Wird die innere 
(mechanische) Dampfung vernachlassigt, wie es meistens bei Schwingungs- 
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rechnungen der Fall ist (um auf der sicheren Seite zu stehen), so wird durch 
die Gleichung A = 0 eine Schwingungsform festgelegt, wie sie gerade im «kri- 
tischen» Zustand des Fliigelflatterns méglich ist. Man sieht aus (21), daB A = 0. 
als Gleichyng zwischen é und 7 (mit den Parametern w und M) durch einen 
Kreis in der (&, 7)-Ebene dargestellt wird. Zur vollstandigen Diskussion setzen — 
wir unter der Voraussetzung kay; + 0: 


A= w havi {(E — Em)® + (1 = tm)? — 77}, (22) 


wodurch die Mittelpunktskoordinaten é,,, 7, und. der Radius 7 des «Grenz- | 
kreises» 2 = 0 definiert sind. Nach Identifikation mit (21) erhalten wir ausge- 
driickt durch die Luftkraftkoeffizienten : 


Rise ake a. 
4 & = ant myt P (23) . 
2 Rays 
kh —k 
be == anr myr F (24) 
i! - Ree 
yee 
a, RR tees (25) 
ayt 


Man sieht aus der Gleichung (22), daB A fiir € = €,, und 7 = 7,, en Extremum 
hat, namlich 
Am = — 0 Ray 7* (26) 


Vor der weiteren Diskussion sei bemerkt, daB beim Ubergang zu einem 
anderen Bezugspunkt nach den Gleichungen (7) die folgenden Transformations- 
gesetze gelten: 


Pe a on 1a Nm >» Ae ae (27) 
und da k,,; nach (7) invariant ist, wird 
PN Sai (27a) 


Es ist also allein €,, keine Invariante, transformiert sich jedoch in anschaulicher 
Ubereinstimmung mit der Gleichung (12a). Die Lage des durch ém festgelegten 
Knotenpunktes relativ zum Profil ist also auch invariant. Bei Beachtung von 
(27) ist es somit gleichgiiltig, von welchem Bezugspunkt ausgehend die Luft- 
kraftkoeffizienten in die Gleichungen (23) bis (26) eingesetzt werden. 

Ein reeller Grenzkreis 4 = 0 ergibt sich nur in dem Falle, wenn der Aus- 
druck unter der Wurzel (25) positiv ist. Wird r imaginar, 7? < 0, so kann fiir 
alle € und 7 nur Dampfung oder nur Anfachung eintreten: die quadratische 
Form A ist definit. Welcher der beiden Falle wirklich vorkommt, sieht man 
am Vorzeichen des Extremwertes 4,,, in (26), wobei also auch das Vorzeichen 
von k,,; eine Rolle spielt. Es sind insgesamt vier Falle zu unterscheiden: 
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I. r? < 0: der Grenzkreis 4 = 0 existiert nicht. 
AL Rays <0, A, < 0: tiberall Dampfung. 
B. Ray; > 9, Am > 0: tiberall Anfachung. 
Wenn der Grenzkreis existiert, so entscheidet /,, das Vorzeichen von A inner- 
halb des Kreises: 
Il. 72> 0: der Grenzkreis 2 = 0 existiert. 
A. Rayi <9, Am > 0: Anfachung innerhalb, Dampfung auBerhalb des 
Grenzkreises. 
B. Rayi > 0, Am << 0: Dampfung innerhalb, Anfachung auBerhalb des 
Grenzkreises. 
Es ist nun eine Besonderheit der Luftkraftkoeffizienten, die spater noch disku- 
tiert werden soll, daB k,,,; immer negativ ist. Darin kommt die bekannte Tat- 
sache zum Ausdruck, daB reine Schlagschwingungen immer gedampft sind, zu- 
‘mindest in Potentialstrémung. Es kommen also nur die Falle IA und ITA in 
Frage: Anfachung tritt nur ein innerhalb des Grenzkreises, wenn dieser existiert ; 
wenn er nicht existiert, so ist die Schwingung immer gedampft, und ein imagi- 
narer Radius r ergibt sich als «hinreichendes» Kriterium fiir Flatterfreiheit: 


potas k t Rian ee r Russ 
ra ( mai 4 ( zn)" eS. (28) 


2 Ray; 2 Ray: ayt 


[Vel. (23), (24) und (25). Wegen ka,; <0 kann die Bedingung (28) nur fiir 
Rmai < 0 erfiillt werden.] 

Wird die innere Dampfung mitberiicksichtigt, so kann man sagen, daB die 
im kritischen Zustand méglichen Werte von & und 7 innerhalb des Grenzkreises 
liegen miissen, denn nur diesen Punkten entspricht eine Schwingungsform, in 
der die Luftkrifte die zur Uberwindung der inneren Dampfung nétigen Arbeit 
leisten k6nnen. 

Die Anwendung der Resultate wird auBerst einfach: die Luftkraftkoeffi- 
zienten sind in Funktion von w und M in verschiedenen Abhandiungen berech- 

- net und tabelliert worden. Mittelpunkt und Radius des Grenzkreises ergeben 

sich nach Einsetzen in die Gleichungen (23) bis (25) in Funktion von w und M. 

Es laBt sich zum Beispiel mit dieser Methode sofort entscheiden, ob fiir ein 
gewisses w und M reine Torsionsschwingungen moglich sind oder nicht. Nach 
dem frither Gesagten sind diese Formen durch 7 = 0 gekennzeichnet. Es kommt 
also nur darauf an, ob der Grenzkreis die Achse 7 = 0 (€-Achse) schneidet. Die 
Bedingung fiir die Existenz von solchen (zwei) Schnittpunkten ist offenbar 


y= ha (29) 
oder nach (25) 
eee 
ee agree (30) 


ayi 
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Wenn die zwei Schnittpunkte existieren, so gibt der zwischenliegende Abschnitt 
der &-Achse Aufschlp8 iiber die méglichen Drehpunkte einer reinen angefach- 
ten Torsionsschwingung. 


Zur Verdeutlichung sei noch eine andere Herleitung der Bedingung (30) ge- — 


geben. Wenn mit einer Riicklage §/ gegen den Bezugspunkt ein Drehpunkt 
teiner Torsionsschwingung mdglich sein soll, k6nnen Aussagen am besten ge- 
macht werden nach Umrechnung der Luftkraftkoeffizienten auf den Dreh- 
punkt als Bezugspunkt. Dies geschieht nach den Formeln (7) mit x = . Ange- 
facht wird eine solche Schwingung dann, wenn das Moment der Luftkrafte, 
welches in Phase mit. der Winkelgeschwindigkeit ist, einen positiven Koeffi- 
zienten hat. Das heiBt aber, daB y,,,; > 0 sein muB, also nach (7) 


Rnai = Rncet + 3 gai te Fm yi) a cs Ra yi BT 


Als Funktion von & ist h,,,; eine Parabel. Es wurde k,,,; < 0 vorausgesetzt, so 


daB fiir = --0o k,,,; <0 wird: reine Schlagschwingungen sind gedampft. — 


Ob diese Parabel iiberhaupt in das positive Gebiet eindringt, kann am besten 
am Vorzeichen seines Scheitelwertes beurteilt werden: es ist 


DR nei 0 
dé 
fiir 
Rows sF Looe Tee 
tee 2k Sm > 


, 


also gerade fiir § = &,,. Der Wert von &,,,,; ist dort 
Rint (€ = Sal a, Rina ie Rayi oF 


Wir haben also die Bedingung fiir Anfachungsméglichkeit einer reinen Tor- 
sionsschwingung : een eon 
was bei negativem k,,,; mit (30) in Ubereinstimmung steht. 

Die bisherigen Anwendungen der Methode, mit Hilfe der Arbeitsberechnung 
die méglichen Flatterformen zu untersuchen, sind schon eingangs erwahnt 
worden. GREIDANUs [2] hat eine sehr vollstandige Diskussion fiir inkompressible 
Strémung gegeben. DuNCAN [16] veranschaulicht die Schwingungsformen mit 
Arbeitsaufnahme und -abgabe durch eine «Flattermaschine». BARTON [3] be- 
handelt Flattern in Uberschallstrémung. In den genannten Arbeiten ist librigens 
das komplexe Amplitudenverhaltnis durch Betrag und Argument gekenn- 
zeichnet: 

Pa lele?. 


Der Zusammenhang mit unserer Bezeichnungsweise ist 


E=|C|cosp, n= |C|siny. 
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Die Beniitzung von § und 7 hat den Vorteil, daB wenigstens 7 eine Invariante 
ist und auch € eine anschaulich verstandliche Abhangigkeit vom Bezugspunkt 
hat. Auch die Kennzeichnung der Schwingungsform durch § und 7 sowie der 
«Grenzkreis» in der (&, 7)-Ebene sind besonders einpragsam. 


4. Quellen fiir die numerischen Werte der Luftkraftkoeffizienten 
a) Unterschallstromung 


Im Sonderfall der Schwingungen in inkompressibler Stromung sind die 
Luftkrafte wohlbekannt; die zu ihrer Berechnung benétigten Funktionen der 
reduzierten Frequenz sind zum Beispiel bei KassNER und Fineapo [1] 
tabelliert. 

Fiir die kompressible Unterschallstromung wurden in den folgenden Rech- 
nungen die nach der Methode von DIETZE [4] aufgestellten Tabellen beniitzt. 
DrevzeE erhalt die Luftkrafte ausgehend vom bekannten Sonderfall M = 0 mit 
Hilfe einer Iterationsmethode, die allerdings um so schlechter konvergiert, 
je naher M an 1 herankommt und je héher die reduzierte Frequenz ist. , 


b) Uberschallstromung 


Fiir M > 1 lassen sich die Luftkrafte in geschlossener Form angeben; es 
sei hier auf die Arbeiten von Possio [5], v. BORBELY [6], ScHwArz [7], HONL [9], 
Tempe und Jann [9], GARRICK und RuBINOW [10], Mires [11] hingewiesen. 
Ausfiihrliche Tabellen finden sich bei TEMPLE und Jann [9], GARRICK und 
RuBrnow [10] sowie bei BARTON [3]; iiber den Zusammenhang der hier verwen- 
deten Luftkraftkoeffizienten mit der Bezeichnungsweise dieser Arbeiten orien- 
tiert die Tabelle 4. Bei den folgenden Rechnungen wurde das von TEMPLE und 
JAHN und von BARTON gegebene Zahlenmaterial beniitzt. 

Es sollen noch die Luftkraftkoeffizienten fiir M > 1 in zwei Sonderfallen 
angefiihrt werden. Zundchst die «quasistationairen» Grenzwerte fiir sehr kleine 
Werte der reduzierten Frequenz w, wobei fiir die Realteile-die nullten, fiir die 
Imaginarteile die ersten Potenzen von @ beibehalten worden sind. Nach den 
Ergebnissen der oben angefiihrten Arbeiten ergibt sich in unserer Bezeich- 
nungsweise mit der Profilmitte als Bezugspunkt: 

4 


Rave =0, Rayd ym? — 1 wo, 
4 2 
k = are TT Ric = eens MO, 
eit = S (ym? — 1)" (31) 
Rar 0 > Rn yi a 0 ) 
M*— 2 
Rema =0 ? Mat 3 (Vm? — 1)° a 


392 Nixo.aus Rort 


ZAMP 


Weiterhin seien die asymptotischen Funktionen angefiihrt, welche die Luftkraft- 


koeffizienten fiir groBe w > oo darstellen. Man findet 


Rayr =O, higyg Pe yp Os 
eee ae 
ae ep Xs 

Rap =O pei ae og 


Ein Vergleich der Werte (31) und (32) zeigt, daB fiir groBe Mach-Zahlen M > oo 
die Luftkraftkoeffizienten sehr angenahert durch die gleichen Funktionen dar- 
gestellt werden sowohl im Grenzfall @ > 0 als auch fiir w + oo. Es liegt die 
Vermutung nahe, da8 dann auch im ganzen Bereich der reduzierten Frequenz 
die einfachen Formeln (32) eine gute Naherung darstellen. Dies ist auch tat- 
sachlich der Fall. Fiir M-Werte nahe an 1 zeigen zwar die Luftkraftkoeffizien- 
ten, die geschlossen durch Besselsche und trigonometrische Funktionen ausge- 
driickt werden kénnen, sehr starke Schwankungen iiber w. Mit gréBer werden- 
der Mach-Zahl glatten sich jedoch die Verlaufe, und von etwa M = 2 an konnen 
die Formeln (32) schon in vielen Fallen als gute Naherung verwendet werden. 
Sie lassen sich nach COLLAR [12] elementar durch eine quasistationare Rechnung 
finden, deren Berechtigung aber — wie TEMPLE und JAHN [9] hervorgehoben 
haben — auf hohe Mach-Zahlen beschrankt ist!). Auch die Schwingungsrechnung 


wird mit den Werten (32) besonders einfach. 


c) Schallstromung (M = 1) 


Im Sonderfall M = 1 sind eigene Rechnungen durchgefiihrt worden*). Nach- 
dem schon in einer kurzen Mitteilung die Methode zur Bestimmung der Druck- 
verteilung bei M = 1 dargestellt worden ist [13], muBten noch die Luftkrafte 
und Momente durch eine Reihe von Integrationen ausgerechnet werden. Diese 
Operationen stehen in vollsténdiger Analogie zu den schon mehrfach dargestell- 
ten Rechnungen fiir M > 1, so daB auf ihre Wiedergabe hier verzichtet werden 
soll. Es seien nur die Resultate angegeben, und es soll gezeigt werden, daB diese 
von den bekannten Ergebnissen fiir M > 1 durch einen Grenziibergang gefun- 


den werden kénnen. 


Die analytischen Ausdriicke fiir die Luftkraftkoeffizienten sind in der Ta- 
belle 1 zusammengestellt. Tabelle 2 gibt die Grenzwerte fiir kleine Werte von a. 


1) In neuester Zeit ist von SAUER [17] eine elementare Herleitung der fiir alle Mach-Zahlen 


giltigen quasistationaren Werte (31) angegeben worden. 


*) Wahrend der Drucklegung dieser Arbeit wurde dem Autor bekannt, daB® eine Abhandlung 


von HEASLET, LoMAx und SpREITER [18] das gleiche Thema behandelt. 


a 


Spe oh Yee 


7 
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Fiir w oo gelten die Formeln (32) mit M= 1. Tabelle 3 enthalt die nume- 
rischen Werte (vgl. auch Fig. 3). 

Beim Grenziibergang gegen M = 1 von der Uberschallseite her geht man 
am besten in zwei Schritten vor. Im ersten — elementaren — Schritt wird 
iiberall M=1 gesetzt, wo dies ohne « Gefahr» méglich ist, das heiBt, wo wegen 


a | 


6 2 


Fig. 3 
Luftkraftkoeffizienten fiir M = 1. Bezugspunkt: Profilmitte. 


- M=1keine GréBe unendlich wird. Als Ergebnis dieses Schrittes erhalten die 
Luftkraftkoeffizienten schon den genau gleichen Aufbau wie in der Tabéile 1, 
nur haben noch die Funktionen G,, Gy, F,, F, eine andere Bedeutung. Aus- 
gehend zum Beispiel von den Resultaten von TEMPLE und JAHN [9] findet man 
sie leicht (unter Beachtung der Zusammenhinge der Tabelle 4) in folgender 
Form: mit der Abktirzung 


qpet a? Ee) 
wird 
ome (est { Jolz) sin (M2) — Jy(2) cos (M2)}, | 
. = (34) 
G, = at { Jo(2) cos (Mz) + J,(2) sin (M2)} , | 
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1 
2) : 
= zu) sin(M zu) du 
= i Jole w) sin (Mu) 


rae qm u) cos (M zu) du 


Fiir M nahe an 1 ist nun z nach (33) groB, und die Bessel-Funktionen konnen — 


durch die Hankelsche asymptotische Darstellung approximiert werden: 


pe V2, cos(e- 2), ley=V sin(e- FZ). G9) 


Der Grenzwert von G, zum Beispiel folgt dann durch Einsetzen von (36) in (34): 


G,= ——— { sin (M z) cos(z _ 4 — cos (M z) sin(z > 3), 


woraus nach einfacher trigonometrischer Umformung folgt 


Z Be I 
| sin |(M—1) 2+ 4], 


und schlieBlich wird nach Einsetzen von z 


Z : wo M qt 
Gyo eee sin(ap 1 =). 
Dies geht fiir M = 1 in den Ausdruck iiber, der in der Tabelle 1 gegeben ist. 
Fiir G, ist die Rechnung ganz analog. 

Beim Grenziibergang fiir F, und F, ist zu beachten, da8 in die Integrale, 
die von der Grenze 0 ausgehen, die asymptotischen Formeln nicht eingesetzt 
werden diirfen. Eine einfache Umgehung dieser Schwierigkeit soll am Beispiel 
von F, gezeigt werden: nach Einfiithrung der neuen Variablen zu = v wird 


E+ le / tay an tat ade 
VM?—1 


I 


seal ge fi Jo(v) sin (M v) dv — / Jo(v) sin (M v) dv 


Das bestimmte Integral zwischen 0 und oo hat den bekannten Wert 1 \VM a 15 
(Vergleiche zum Beispiel die Untersuchungen iiber die Integrale vom Typus 


(35) von ScHWARz[14].) Im zweiten Integral kann die asymptotische Formel 
(36) eingesetzt werden: 


R=y [VS  cos(v — ©) sin (Mv) dv. 


(35) 
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_ Nach einer einfachen Umformung der trigonometrischen Funktionen wird 


Pariser pS 


ope YEP sin [we — 1 = dv 
— VT fin lar + ee alae ae 


Mit den neuen Variablen (M — 1) v =¢ und (M+ 1) v=s erhalt man 


Cc 
dt yi ; m ds 
ees 


14 
7) 2 nt M hae ee 
V2at ee V2as 


Fiir M = 1 wird nach Einsetzen von z in Ubereinstimmung mit dem Ausdruck 
in der Tabelle 1: 

ire) dt w/2 
Ped -| (sin¢ + cost) —=—— =f (sin ¢ + cos?) et ay 
y2 ut " y2 at 


w/2 


- denn es ist der bekannte Grenzwert der Fresnelschen Integrale zwischen den 
Grenzen 0 und oo: 


co co 
: at dt 1 
sin ¢ ———._ = cost ———_- = =. 
V2at V2xt - 
0 0 


Die Luftkraftkoeffizienten fiir M = 1 sind in Fig. 3 aufgetragen; sie zeigen 
iiber w einen glatten Verlauf. Die benachbarten Verléufe M > 1 hingegen 
weisen iiber w rasche Schwankungen auf. Dazu sei bemerkt, da beim Grenz- 
iibergang gegen M = 1 von der Uberschallseite her Glieder vom Typus 


ae 


(M — 1) sin 7 
zum Verschwinden kommen. Der Durchgang bei M= 1 erfolgt also tiber M 
mit immer rascher und kleiner werdenden Schwankungen. Fiir groBe w-Werte 
ist allerdings das Schwankungsgebiet praktisch sehr klein. Bei Verkleinerung 
yon w wird man aber zu erhodhten Schwankungen und schlieBlich im Grenzfall 
@=0 zu Diskontinuitaten der Luftkraftkoeffizienten bei M=1 gefiihrt. 
Mathematisch sind diese wegen Divergenz der asymptotischen Formeln bei 
@ = 0 verstandlich. Es kommen damit bekannte Anderungen der Stromung 
beim Schalldurchgang zum Ausdruck: zum Beispiel ist der stationare Kraft- 
angriffspunkt im Unterschallbereich der vordere Viertelpunkt, bei Uberschall 
die Profilmitte, wahrend fiir M= 1 die (unendliche) stationare Luftkraft in 
ebener Strémung nach linearisierter Theorie in 1/3-Profiltiefe angreift [13]. 


a yf Ng tage 
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Wenn man aber das Gebiet in der Nahe von w = 0 von der Beniitzung aus- © 


schlieBt (wozu noch weitere Griinde angefiihrt werden sollen), so kénnen die 
Luftkraftkoeffizienten fiir M = 1 gewissermaBen als Mittelwerte im schallnahen 
Gebiet angesehen werden. Der erforderliche Rechenaufwand fiir eine Schwin- 
gungsrechnung in der Umgebung von M=1 ware dadurch erheblich ver- 
mindert. 

Bei Verwendung der Luftkraftkoeffizienten fiir M = 1 sind vor allem zwei 
physikalisch begriindete Einwande zu beachten: 

1. Es kann bezweifelt werden, daB die Ergebnisse einer linearisierten 
Theorie fiir = 1 brauchbar sind. Die Untersuchungen von Lin, REISSNER 
und TsIEN [15] haben gezeigt, daB die Linearisierung selbst fiir M = 1 zulassig 
wird, wenn in der Potentialgleichung die linearen nichtstationdren Glieder 
gegen die vernachlassigten quadratischen Glieder tiberwiegen, das heiBt, wenn 
die Schwingung gentigend rasch, die reduzierte Frequenz geniigend hoch ist. 
Es kommt darauf an, da8 der Wert w/t?/? groB genug wird, wobei t ein « Dicken- 
parameter» ist (prozentuelle Dicke oder Anstellwinkel des Profils). 

2. In der naheren Umgebung von M = 1 versagt die sogenannte Streifen- 
methode, das hei8t wenn der wirkliche Fliigel endliche Spannweite hat, so ist 
die Betrachtung der ebenen Strémung eine unzulassige Idealisierung. Die 
Ubertragung der folgenden Resultate auf den dreidimensionalen Fall ist in der 
Umgebung von M = 1 nicht ohne weiteres zulassig. 


Tabelle 7 
Die Luftkraftkoeffizienten fiir M = 1. Bezugspunkt: Profilmitte 
Ray, = —4 @ (G, — B) 
Rays =a (Gy aie F,) 
Et 
Paar = 4(G, +R) +2(2=G,] 
bagi = —4 (G— FR) + 2(G, — 4) 
cai Fy 
hye = ~2 (2 — Ga} 
ene ys 
Bye ‘ 2 ( 1 “| 
a) 2 va 4 
bnar = — 3 (GA-F)-S(G-S)+2G 
o Deoidise 4 
bn = 3 (+A) +o(2-G)-+G, 
dabei ist 
1 wo : 1 
(é 1 i (00) : (20) 
1 eRe (cos 7 + sin a G,= == (cos ees al 3) 
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@/2 F @/2 
: t ; 
B as (cosé + sin?) = Hira (cos¢ — sin #) Lie 
3 V at ; y2 ut 
Tabelle 2 


Naherungsformeln zu den Ausdriicken der Tabelle 1 fiir kleine Werte der redu- 
zierten Frequenz w 


a 4 a as — 48 

Veh = et glo Vakg =— at glo 

sac hae ae 1 = — 28 2s 

Vit hing, = — Vo —-= oVo Vi ky: = — = Vo+~oVo 

= LOE EGU | oh a eo 4 4/- 
Tabelle 3 


Die numerischen Werte der Luftkraftkoeffizienten fir M= 1; 
Bezugspunkt: Profilmitte 


@ —Rayr —Rayi Raar —Fagi 
0,04 0,4424 0,4604 11,5863 10,9846 
0,08 0,6129 0,6640 8,4090 7,5580 
0,12 037353 0,8291 7,0443 6,0023 
0,16 0,8315 0,9759 6,2566 5,0536 
0,2 0,9100 sta ii lp: SoTL 4,3924 
0,4 1 ism Eve 1,7219 4,5690 2,6719 
0,6 1,2514 2,2971 4,1662 1,8519 
0,8 1,2689 2,8751 3,9985 1,3407 
iL 1,2309 3,4690 3,9355 0,9849 
1,2 1,1529 4,0843 3,9278 0,7234 
1,4 1,0455 4,7236 3,9521 0,5258 
1,6 0,9164 5,3880 3,9954 0,3752 
1,8 0,7718 6,0779 4,0500 0,2607 
eZ 0,6168 6,7929 4,1109 OL 752 

Tabelle 4 


A. Zusammenhang mit der Bezeichnungsweise von TEMPLE und JAHN [9]. 


Bezugspunkt: Profilmitte. 
Rage = 24, 
k 


ag 


aA 


—2a]; ae ole 


—Rmyr —Rm yi Rmar —Fmai 
0,0761 0,0743 1,7905 1,9710 
0,1089 0,1038 1,2026 1,4579 
0,1349 O12595 0,9303 1,2429 
0,1575 0,1430 0,7611 12221 
0,1779 0,1578 0,6411 1,0446 
0,2643 0,2074 0,3149 0,8848 
0,3378 OFZ537 0,1469 0,8433 
0,4045 0,2454 0,0346 0,8361 
0,4665 0,2460 —0,0491 0,8434 
0,5243 0,2374 —0,1149 0,8580 
0,5781 0,2209 —0,1682 0,8765 
0,6278 0,1976 +0,2121 0,8976 
0,6734 0,1681 —0,2484 0,9202 
0,7144 0,1332 -—0,2782 0,9439 
=—2m, gar = 4 My 
=-—20m; Rei = 2 O Mg, 
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B. Zusammenhang mit der Bezeichnungsweise von BARTON [3]. i 

Bezugspunkt: Profilvorderkante. ; 

nye = AO Py Ria ZOE, Kiggg = 20" Mg Rap =-2 0" Msg $ 
yi = 20° RO k= 207k kimyi = 20° Ms mai = 20" My, 


C. Zusammenhang mit der Bezeichnungsweise von GARRICK und RuBINOW [10]. 
Bezugspunkt: Profilvorderkante. 


ae 
/ ie f es 2 y, / —— 2 , , ast vy 
Raye = —4 ol, k,,=o'L, (eS = M, Ruae = 3 M3 
2 oe ae ne Ta ees 
Rayi = —2.07 L, Rai NILES) myi © 2 Mot hee 


Bemerkung: w ist die mit der vollen Profiltiefe errechnete reduzierte Fre- 
quenz. 
5. Resultate 


Mit Hilfe des genannten Quellenmaterials wurden die Mittelpunktskoordi- 
naten €,,, mm und der Radius 7 des Grenzkreises nach den Formeln (23) bis (25) 


ii 


08 


06 


0 04 OP 12 ow 46 — 0 O4 08 42 16 w 2 
Fig. 4 

Mittelpunktskoordinate €,, des Grenzkreises tiber der reduzierten Frequenz qm fiir verschiedene 

Mach-Zahlen M. 


ausgerechnet. Die Ergebnisse sind in den Fig. 4 bis 6 zusammengestellt. Auf- 
getragen wurden €,, und die (mit m multiplizierten) Werte w 7,,,und w r tiber 
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w fiir verschiedene Mach-Zahlen. Berechnet und in Fig.7 dargestellt ist weiter- 

hin die GroBe Z,, nach (26), welche ein MaB fiir die gré8tmégliche Anfachung 
- gibt. : 

Zuniachst sei das auffallende Resultat hervorgehoben, daB (mit Ausnahme 

des inkompressiblen Falles M = 0) der Radius 7 und damit auch /,, bei einem 


16 wie 
Fig. 5 
Mittelpunktskoordinate 7,, des Grenzkreises iiber w fiir verschiedene M. 


|" 


Ob 


95 9 04 08 42 


05 


1} G4 Lis THF Palaad lb hy 08 ri eee (Maas 


Fig. 6 
Radius r des Grenzkreises iiber w fiir verschiedene M. 


endlichen, von der Mach-Zahl abhangenden Wert w = * zu Null wird; fiir 
groBere Werte von w wird 7 imaginir, /,, negativ. Es tritt also fiir w > w* der 
vorhin besprochene Fall ein, daB das hinreichende Kriterium (28) zur Flatter- 


oe 
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freiheit erfiillt ist. Den Verlauf von w* iiber M zeigt Fig. 8. Fiir die Flatter-_ 
rechnung ergibt sich die wichtige praktische Konsequenz, daB die kritischen q 
GréBen nur unterhalb der Grenze w* gesucht werden miissen. Auch bei der 
Aufstellung von Tabellen der Luftkraftkoeffizienten kann man sich vom Stand-_ 
punkt des Flatterproblems aus mit zwei Freiheitsgraden auf das Gebiet w <ar" 


hy aT er a eee 


Fig. 7 
Maximaler Leistungsbeiwert A,,,. 


0 O4 08 42 Lo Gr Ke 


beschranken. Fiir M = 1 wurde zum Beispiel w* = 1,809 gefunden. — Wird die 
Zahl der Freiheitsgrade durch Einbeziehung der Ruderschwingungen erhoht, 
so verschiebt sich allerdings die Grenze weiter nach oben. 

Allein im Sonderfall M = 0 findet man w* = oo; Flattern ist also prinzipiell 
bis zu beliebig hohen reduzierten Frequenzen méglich. Die fiir M = 0 berech- 
neten Luftkraftkoeffizienten fallen aber gewissermaBen aus der Reihe, da bei 
ihrer Berechnung nicht nur die Geschwindigkeit U klein gegen die Schallge- 
schwindigkeit vorausgesetzt war, sondern auch die Schallgeschwindigkeit = co 
gesetzt wurde. Dies kann bei hohen Frequenzen (das heiBt groBen w-Werten) 
selbst bei kleinem U zu unrichtigen Ergebnissen fiihren. Es ist wahrscheinlich, 
da bei Beriicksichtigung der endlichen Schallgeschwindigkeit auch im Grenz- 
fall U > 0 ein endliches w* gefunden wird. Fiir M = 0,5 und M = 0,6 gestatten 
es die vorhandenen Werte von DIETZE nicht, auf w* zu extrapolieren. Fiir 
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M = 0,7 findet man w* = 1,78; der Wert wurde als Nullstelle von A,, inter- 
poliert zwischen w = 1,4 und w = 2 (Fig. 7), so daB die letzte Stelle unsicher 
ist. Es sei auch daran erinnert, daB die Rechnung nach Dietze bei hohen 
- M- und w-Werten schlecht konvergiert. 


od ae 


G5 


a 05 1 (AMO AIZ 
Fig. 8 


Verlauf der Grenze w* fiir Flattern mit zwei Freiheitsgraden tiber der Mach-Zahl. w%. = Grenze 
fiir Flattern mit einem Freiheitsgrad (Torsion), 


Im Uberschallgebiet wurden mit den Tabellen von TEMPLE und JAHN die 
folgenden Werte gefunden: 


Interpoliert zwischen 


M w* w< w* o> w* 
2 1,69 (1,53 und 1,83) 
1,4 1,42 (0,98 und 1,47) 
1,6 1,43 (1,22 und 1,82) 
ats) 15d (1,38 und 2,07) 
io) 155 (155\ Sund 1,65) 


Im Grenzfall M oo lassen sich durch Einsetzen der Werte (32) in die Formeln 
(23) bis (26) leicht die folgenden analytischen Ausdriicke finden: 


m = O (Profilmitte) , 
i 1 1 1 1 w? 
2 e-f. a -0-9) 


Aus den Formeln (37) folgt @* = V3 = 1,732 fiir M >. 
Geometrisch besagt unser Ergebnis, da8 Flattern nur eintreten kann, wenn 


nach der Beziehung (4a) 


(37) 


ZAMP 1/26 
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ist; zum Beispiel fiir w* = V3 muB die Wellenlainge L gréBer als 3,62 Fliigel- 
tiefen J sein. Dies ist in gewissem MaBe auch anschaulich zu verstehen: fiir 
Wellenlangen, die zu klein sind verglichen mit der Fliigeltiefe, kann sich das 
Profil nicht mehr so «durchschlangeln», daB die Schwingung ungedampft oder 
' gar angefacht wird (vergleiche Fig. 2). 

Als nachstes seien die quasistationaren Grenzwerte fiir kleine w diskntees 
In stationarer Stroémung w = 0 ist fiir M< 1 


LGC 
aie | 
und fir M >1 (38) 

ea rae | 


Mit den entsprechenden Kraftangriffspunkten als Bezugspunkt (M <1: 
1/4-Punkt, M > 1: Profilmitte) findet man durch Vergleich mit (3) fiir m = 0: 


Ba Moods. k ae 


Mg be eee Sah Pee 
aar yi- We ~M?—1 


(39) 


alle anderen Luftkraftkoeffizienten = 0. 

Nun denken wir uns «= const = %, daneben soll aber das Profil eine 
«langsame» Schlagschwingung ausfiihren: y = y, e'”’. Quasistationar kénnen 
dann die Luftkrafte bekanntlich durch Beriicksichtigung eines scheinbaren 
Anstellwinkels «; = —y/U in den stationaéren Formeln (38) erfaBt werden: 
a ist mit « + «> zu ersetzen. Es ist 


to, = aaa a ; 
also die entsprechende Amplitude 
as, = — V9 (40) 
und nach (38), (39) und (40) 
Ca, = Pans (x —141@ Ae), (41) 
nach Vergleich mit (3) folgt 
Rays = — O Ry ae (42) 


als quasistationar allgemein giiltige Grenzbeziehung. Zunachst sieht man, daB 
Ray: nach (42) negativ ist. Es ist anschaulich leicht einzusehen, daB quasista- 
tionar die reine Schlagschwingung immer gedimpft wird, was gleichbedeutend 


wt = 
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ist mit negativem k,,;. Das groBe Zahlenmaterial mit exakten Werten fiir kay; 
zeigt fiir alle m und M, daB k,,; < 0, auch genau gerechnet, immer zutrifft. 
Die Beziehung (42) gestattet schon, den fiir alle Mach-Zahlen giiltigen und 


- fiir die Schwingungsform entscheidenden Grenzwert von 7, fiir w—> 0 nach 
(24) anzugeben. Der Grenzwert von k,,,, wird mit dem stationaren Kraft- 


angriffspunkt als Bezugspunkt gleich Null, und man erhalt nach (24) und (42): 


; il : 1 
ie ees AD LR Mei et (43) 

Fiir den Grenzwert von &,, braucht man nach (23) noch die quasistationaren 
Werte von k,,,,; und ky,;. Aus (41) folgt kya: = Rmyi = 0 fiir w > 0; da aber 
auch k,,,; gegen Null strebt fiir m > 0, ist &,, mit den bisherigen Ergebnissen 
noch unbestimmt. Fiir M< 1 seien nur die Ergebnisse einer weitergehenden 


Untersuchung angefiihrt: es wird 


limé,, = —co, 


o=0 


und zwar logarithmisch unendlich mit . Fiir M = 1 folgt aus den Werten der 
Tabelle 2 


- 1 
o=0 y 
wahrend sich fiir M > 1 aus den Formeln (31) ergibt: 


1 
lim em = —Z Ue — 1) 


w=0 


mit der Profilmitte als Bezugspunkt. 
Fiir y findet man (mit der bemerkenswerten Ausnahme M = 1) , daB fiir 
@ > 0 im Ausdruck (25) die unendlich werdende GroBe 7), gegen alle anderen 
tiberwiegt, so daB der Grenzwert gilt (M + 1): 
: 1 
im? = tim) 7,,| =>—- (44) 


2 
wo=0 wo=0 


Im Sonderfall M = 1 ergibt sich quasistationar fiir 7 eine Diskontinuitat : nach 
den Werten der Tabelle 2 wird 


iy geo CE (45) 
V2 w=0 


o=0 
Man sieht, daB 7,, und r fiir alle Mach-Zahlen, aber auch €,, tiir M <1 gegen 


Unendlich streben mit @ > 0. Aus diesem Grunde sind in den Fig. 5 und 6 
die endlich bleibenden Werte w7,, und wr aufgetragen. Es ist auch fiir sehr 


egos 
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kleine w-Werte giinstiger, den Grenzkreis in einer (w €, w y)-Ebene zu unter- 

suchen. Geometrisch bedeutet die Multiplikation mit wm die Einfiihrung der — 

Wellenlange L statt / als Bezugslinge: nach (4a) und (8) ist 
Vo Vo 


i 
oC=2n7- oe = 290 ree a 


Auch die Betrachtung einer «Grenzellipse» in einer (§, w7)-Ebene kann 
praktisch vorteilhaft sein; fiir die Schwingungsform ist nach (16) direkt die 
GréBe w 7 maBgebend. 


Fig. 9 
Grenzkreise fur jm = 1 und verschiedene Mach-Zahlen. 


Bei der weiteren Verfolgung der Verlaufe von &,,, 17, und 7 tiber w zeigt sich, 
daB besonders die Werte von w7,, und wy fiir verschiedene Mach-Zahlen 
nirgends um gréBere Betrage voneinander abweichen (Fig. 5 und 6). Die Ver- 
laute fiir Uberschall-Mach-Zahlen nahe an 1 zeigen rasche Schwankungen, wie 
es nach dem entsprechenden Verhalten der Luftkraftkoeffizienten zu erwarten 
war; die Abweichungen von den Werten fiir M = 1 erreichen aber nur relativ 
kleinere Betrage. Man darf also schlieBen, da8 im allgemeinen die kritischen 
Schwingungsformen nur wenig von der Mach-Zahl abhangen. Als Beispiel sind 
in Fig. 9 die Grenzkreise fiir w= 1 und fiir die Mach-Zahlen M = 0; 0,7; 1; 1,4 
eingezeichnet. Lage und Gr6éBe der Kreise sind nicht stark verschieden. Es ist 
von besonderer Bedeutung, daf alle diese Kreise ganz in der Halbebene y < 0 
legen. Unter dieser Bedingung haben namlich die gleichen konstruktiven Richt- 
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linien zur Verhiitung des Flatterns ihre Giiltigkeit, die schon bei der Diskussion 

des inkompressiblen Falles aufgestellt wurden. Die MaBnahmen lauten: 

_ 1. Biege- und Torsionsfrequenz der Standschwingung sollen méglichst verschie- 
den sein. Da bei den gefahrlichen Flatterformen beide Freiheitsgrade mit- 
spielen, ist es vorteilhaft, wenn die beiden Frequenzen héchstens durch 
groBe Luftkrafte zueinander gestimmt werden kénnen. 

2. Der Schwerpunkt soll méglichst vorne liegen. Die Schwingungen mit 7 < 0 
sind namlich dadurch gekennzeichnet, daB zum Beispiel beim Herunter- 
schwingen die Profilhinterkante gegentiber der Vorderkante zuriickbleibt, 
wie aus Fig. 2a ersichtlich ist (im Gegensatz zu Fig. 2c). Eine derartige 
gefahrliche Schwingungsform wird durch einen hinten liegenden Schwer- 
punkt begiinstigt, so da8 umgekehrt die Schwerpunktsvorlagerung die 
Flattergefahr vermindern wird. 

3. Die beste Lage der elastischen Achse wird bekanntlich durch die statische 
Bedingung des Auskippens bestimmt und ist méglichst vorne am giinstig- 
sten. Vom Flatterstandpunkt aus wird allerdings die gefahrliche Form 
(Fig. 2a) durch vorne liegende elastische Achse begiinstigt; weit hinten 
liegende elastische Achse bedeutet hingegen Schwachung der Torsionsfeder, 
so daB man hier auf diesem qualitativen Wege zu keiner weiteren Empfeh- 
lung gelangen kann. 

Es ist auch anschaulich einzusehen, daB die in Fig. 2a dargestellte Schwin- 
gungsform gefahrlich ist. Die Luftkraft wird namlich naéherungsweise (quasi- 
stationar) bestimmt durch den Anstellwinkel relativ zur absoluten Bahn des 
Profils. Bei der Form nach Fig. 2a ergibt sich dadurch eine Kraft nach unten 
beim Herunterschwingen und nach oben beim Heraufschwingen, also Anfa- 
chung. 

Die groBtmégliche Anfachung ergibt sich fiir die Schwingungsform, die dem 
Mittelpunkt des Grenzkreises entspricht. Man versteht, daB die Mittelpunkte 
fiir alle M und o in der negativen Halbebene 7 < 0 liegen. Nach Fig. 5 wird 
sogar die Abweichung vom quasistationdren Grenzwert (43): @ %m = —9,5 
nicht sehr groB fiir alle M und o. 

Die Ergebnisse der Flatterrechnung von GARRICK und RuBinow [10] fiir das 
Uberschallgebiet bestatigen die bisherigen qualitativen Folgerungen. Betref- 
fend der Schwerpunktlage finden GARRICK und RuBINow, daB es eine Grenz- 
lage gibt, bei der durch weiteres Vorriicken des Schwerpunktes das Flattern 
iiberhaupt unméglich wird. Die Grenze liegt etwa bei der Profilmitte fiir solche 
Torsionsfrequenzen, die viel hoher sind als die Biegefrequenz; bei Annaherung 
der beiden ist’ die Grenze eher durch die elastische Achse gegeben. Wenn der 
Schwerpunkt hinter diese Grenze gelangt, wird aber die Flattergefahr meist 
sehr plétzlich erhéht. 

In einem Falle gelangen jedoch Garrick und RuBinow zu keiner eindeu- 
tigen Flattergrenze durch Schwerpunktvorlagerung, namlich bei der kleinsten 


} 
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untersuchten Mach-Zahl (M = 10/9), die M = 1 am nachsten liegt. Der Mach- z 
Bereich nahe an M = 1 ist noch durch eine andere Erscheinung gekennzeichnet, ~ 
die von verschiedenen Autoren schon hervorgehoben wurde: es besteht die 
Moglichkeit von Flattern mit nur einem Freiheitsgrad, das heiBt in der Form — 
reiner Torsionsschwingung. Nach unserer Betrachtungsweise hangen diese bei- — 
den Erscheinungen zusammen. Wie schon friiher gesagt wurde, kénnen reine ~ 
Torsionsschwingungen angefacht werden, wenn die §-Achse vom Grenzkreis 
geschnitten wird. Ein Teil des Grenzkreises wird dann in die positive Halb- 
ebene 7 > 0 hereinragen; die gefahrlichen Schwingungsformen, welche diesen 
Punkten entsprechen, werden durch die MaBnahme der Schwerpunktsvorlage- 
rung nicht verhindert, ja sogar begiinstigt. Die Frage des Flatterns mit nur 
einem Freiheitsgrad erfordert also besondere Beachtung. 

Es zeigt sich, daB die zum Schnitt mit der é-Achse notwendige Bedingung 
(29) r > || nur in einem beschrankten Bereich von w und M erfiillt ist. Im 
quasistationaéren Grenzfall sind mit Ausnahme von M= 1 die Grenzwerte r 
und |7,,,| nach (44) einander gleich, das heiBt, die &-Achse wird gerade beriihrt. 
Bei M = 1 gibt es eine Diskontinuitat: nach (45) ist fir w >0:7r> |ym\, 
die é-Achse wird geschnitten. Mit wachsendem w wird aber der UberschuB von 
rv gegen || geringer, und bald wird eine Grenze erreicht — sie soll w} genannt 
werden —, wo der Grenzkreis die €-Achse nur noch beriihrt. Fir M@ = 1 wird 
r= |%m| fiir wo} = 0,465. Bei gréBeren reduzierten Frequenzen ist reine Tor- 
sionsschwingung nicht mehr moglich. 

Im Gebiet M <.1 wachst von w = 0 aus gehend 7 zunachst etwas starker 
als |7%m|, doch bleibt der UberschuB sehr gering und hort bei einer niederen 
Grenze w# auf. Es ergeben sich fiir w* die Werte: 


iy 
q 


| 


M ox 
0) 0,088 
0,5 0,15 
0,6 0,18 
0,7 0,21 


Bei Annaherung an M = 1 wird also w* gréBer; auch von der Uberschall- 
seite her zeigt es sich, daB reine Torsion und Formen mit 7 > 0 besonders in 
der Umgebung von M = 1 auftreten kénnen. Die Werte von w*% sind fiir M > 1: 


M OP 
nial 0,50 
ibe 0,66 
125 0,69 
Ales} 0,68 
1,4 0,62 


Der Verlauf von w# iiber M ist in Fig. 8 eingetragen. Uber den Durchgang 
dieser Kurve bei M = 1 wissen wir, daB er mindestens von der Uberschallseite 
her durch immer raschere und kleinere Schwankungen erfolgt. — Um auch ein 
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Bild zu erhalten von der GréBe des Grenzkreisgebietes, der in die positive 
Halbebene 7 > 0 zu liegen kommt, ist in Fig. 10 der Verlauf von o (7 — | jm]) 
_ uber w fiir verschiedene Mach-Zahlen aufgetragen. Mit Ausnahme der naheren 
- Umgebung von M=1 (besonders tiberschall- 
seitig) sind die Schnitte sehr schleifend. 

In Fig. 8 fallt die w%-Kurve von etwa 
M = 1,4 an steil herunter und erreicht den 
Wert 0 bei der Mach-Zahl M = 5/2 = 1,581. 
Dieser bemerkenswerte Grenzwert laBt sich 
leicht herleiten, wenn man die Grenze fiir reine 
Torsionsschwingung nach der Bedingung (30) 
mit den quasistationaren Luftkraftkoeffizien-  y.uaut von ea an abe Oe 
ten (31) sucht. Er ist auch in den Arbeiten verschiedene M. 
von GARRICK und Rusinow [10] und von 
BarTON [3] enthalten. Fiir gréBere Mach-Zahlen kénnen iiberhaupt keine rei- 
nen Torsionsschwingungen auftreten. 

Die relativ niedrige Grenze der reduzierten Frequenz fiir reine Torsion 1aBt 
schlieBen, daB man dieser Art der Flattergefahr durch geniigend hohe Torsions- 
steifigkeit begegnen kann. Betrachten wir zum Beispiel eine Anstrom-Mach- 
Zahl M = 1,25 (also etwa U = 375 m/s), fiir den sich der héchste Wert w7 = 0,69 
ergeben hat. Nach (4) ist y= @ U/l, bei einer Profiltiefe / = 2m entspricht 
somit dem Wert von w* eine Kreisfrequenz v¥ = 130 rad/s oder die Frequenz 
f* = 20,6 Hz. Wenn bei geniigender Torsionssteifigkeit diese Frequenz (auch 
unter Einwirkung der Luftkrafte) nicht unterschritten werden kann, ist das 
ganze Gebiet 7 > 0 aus dem Gefahrenbereich ausgeschlossen. Die vollstandige 
Diskussion dieses Problems wird erschwert durch den Umstand, da8 fiir M 
nahe an 1 und kleine w-Werte die frither erwahnten physikalischen Einschran- 
kungen gelten. 

Im Rahmen dieser Theorie, welche ebene Strémung und Zulassigkeit der 
Linearisierung voraussetzt, kann man jedoch mit Unterstiitzung der quanti- 
tativen Ergebnisse von GARRICK und RuBINOW zu einer einfachen Empfehlung 
gelangen, wie ein Fliigel am besten fiir alle Mach-Zahlen flattersicher gebaut 
werden soll: durch passende Vorlagerung des Schwerpunktes wird die Flatter- 
gefahr auf die Umgebung von M = 1 beschrankt; diese Gefahr ist fiir die n6tige 
Torsionssteifigkeit maBgebend. 


Anhang: Beziehungen zur volistandigen Lésung des Flatterproblems 


Im folgenden sei noch untersucht, was die bisherigen Uberlegungen vom 
algebraischen Standpunkt fiir die eingangs beschriebene vollstandige Aufgabe 
der Frequenzbestimmung bedeuten; von den Einzelheiten des mechanischen 


| 
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Problems wird abgesehen. Die Aufgabe lautet bei  Freiheitsgraden: es a 
die Eigenwerte A ( = 1/v?) des Gleichungssystems 


Sipe aaee Gf 21,277 2p a) Soe 
k=1 


bei komplexem c;, anzugeben, und es sind Bedingungen fiir die c;;, zu finden, 
unter denen A rein reell wird. Die Gleichungen (46) sind etwas verallgemeinert — 
gegeniiber dem iiblicherweise betrachteten Fall durch Beriicksichtigung der 
Koeffizienten g;, diein unserem Problem wesentlich reell und positiv sein missen. 

Wir beschranken uns im weiteren wieder auf m = 2; die Verallgemeinerungs- 
méglichkeit wird leicht ersichtlich sein. Alle GréBen seien in Real- und Imagi- 
narteil zerlegt : 


Cin = Ayr tt dyn, j= Xi +1Y;, A= 6 +16: (47) 
(47) in (46) eingesetzt, gibt fiir m = 2 vier reelle Gleichungen: 


yy %1 + Ayg %q — Oy Vr — Oye Vo = 0 81% — OHM» 
Ay, %1 + Ogg X_ — Day Vy — Oop Vo = O Bo Xq — O Ba Vo, 


(48) 
By % + Oyy %y + Ay Va + Aye Ve = OU t+OR %, 


Boy %y + Dyn %y + Ay, Vy + gy Vo = O Bo Vo + O Bo Xo - 


Gesucht sind die Bedingungen fiir die a;; und 6;,, unter denen A rein reell, das 
heiBt o = 0 sein kann. Um daritiber zu entscheiden, multiplizieren wir nach 
klassischem Vorbild die Gleichungen (48) der Reihe nach mit — y,, — vo, %1, %» 
und addieren sie. Man kann sich davon tiberzeugen, daB diese formale Operation 
der Berechnung der Arbeit entspricht. Es wird so @ eliminiert, und man findet 
nach leichter Umformung: 


Dar (44 + Yt) + (Oya + Poy) (4 Xe + Va V2) + boa (23 + V3) | 


(49) 
+ (4, — a1) (%1 Yo — Xe Ya) = 0 (81 *7 + Bi VE + Bo XB + Bo V3) - | 


Wenn nun die g; positiv sind, so kann der quadratische Ausdruck auf der 
rechten Seite von (49) auBer dem trivialen Fall x; = y; = 0 nicht verschwinden 
und ist immer positiv. Das Vorzeichen von o wird durch die linke Seite von 
(49) bestimmt. Ist nun 


Ayy = My, Oy = —dby, 0 = dy =0, 
allgemeiner : 
Ayn = Anz, 0 


a — dp; , (50) 


so verschwindet die linke Seite identisch, und es wird immer o = 0. Bei Er- 
fillung der Bedingung (50) ist c,, eime sogenannte Hermitesche Form, fiir 


ik 
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_ welche wir den klassischen Beweis wiederholt haben, daB die Eigenwerte A 
immer reell sind. 
Ist nun die Form c;, nicht hermitesch — und dies ist beim Flatterproblem 
der Fall —, so muB die auf der linken Seite von (49) stehende quadratische 
Form indefinit sein, das heiBt, fiir gewisse Werte von x; und y, (auch auBer 
_%; = y; = 0) verschwinden kénnen, damit o = 0 sein kann. Die Koeffizienten 
dieses quadratischen Ausdruckes enthalten nur 


Ain—@yj und 05,4 04; 
(somit auch die 0;,). 

Beim mechanischen Problem ohne Luftkrafte und ohne Dampfung ist die 
Form c;;, sogar reell symmetrisch, das heiBt a; = a;,; und b;, = 0. Es ist ver- 
standlich, daB dann beim Flatterproblem die linke Seite von (49) nur die Luft- 
kraftkoeffizienten enthalten wird. Die a;; treten in (49) tiberhaupt nicht auf; 
dem entspricht die Tatsache, daB in der fritheren Gleichung (21) Ray, und 
Rmar Dicht vorkommt. 

Fiir n = 2 wird die Diskussion der Gleichung (49) sehr einfach. Wegen der 
Homogenitat des Problems darf es nicht auf die Werte von z, und z, selbst an- 
kommen, sondern nur auf deren Verhaltnis ¢: 


oat Wes 
Ete ee Seb Ait 
oder anders ausgedriickt: mit der Normierung z, = 1, das heiBt 


hye l,. Vg oe (51a) 
wird z, = ¢, das heiBt 
Hae, A= 1. pth) 


Mit (51a) und (51b) wird die Gleichung (49): 
6, (67 + 74°) + (x + bey) € — (@yq — An) 9 + Fog = 7 (g, £7 + #1. + 8)- (52) 


o = 0 ist in der (&, 7)-Ebene die Gleichung eines Kreises; ist dieser imaginar, 
so kann o nicht verschwinden. Bei reellem Kreis kann o = 0 sein, doch ist die 
Realitat des Kreises o = 0 nur eine notwendige, aber keine hinreichende Be- 
dingung dafiir, daB der Eigenwert A von (46) reell ist. 
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Summary 


On the basis of energy considerations a survey is given of the possible forms 


of oscillations for flutter with two degrees of freedom in a plane compressible 
potential flow. It is found that the possible forms do not much depend on Mach 
Number in the whole range from 0 to oo, with some exceptions in the neighbour- 
hood of M = 1, where forms with only one degree of freedom (pure torsion) may 
occur. Conclusions are drawn for methods of preventing flutter. Limits of the 
reduced frequency (depending on M) are given for the possibility of flutter in 


two 


and one degree of freedom. 


Special care was given to the case M = 1, for which analytical expressions and 


numerical values of the derivatives are presented. 


(Eingegangen: 27. 3. 1950.) 
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Kurze Mitteilungen — Brief Reports - Communications bréves 


Zur Multiplikation reeller Zahlen 


Von WALTHER EBERL, Wien’) 


Im folgenden wird das Verfahren der abgekiirzten Multiplikation so erweitert, 
wie es fiir die Zwecke des Rechnens mit Rechenmaschinen erforderlich ist. Fur 
diese Art der abgekiirzten Multiplikation werden Genauigkeitsschranken abge- 
leitet, die diejenigen der gewohnlichen abgekiirzten Multiplikation als Spezialfall 
enthalten. 

Dem Vernehmen nach ist diese Methode den Fachleuten der USA., die mit 
elektronischen Rechenmaschinen arbeiten, bereits bekannt; doch da Publika- 
tionen dariiber nicht bestehen, erscheint die folgende Notiz nicht iiberflissig. 

Hat man zwei reelle Zahlen 

=: Fe —co 

a= Ya,10! und b= Jb, 10 

m n 
zu multiplizieren, so geschieht das am besten nach Art der Multiplikation von 
Reihen: 

=O) : wn 

ab= ¥c,10', wo ¢,= 3 4,),_; 

m+n j=m 
ist. Wahrend die a; und 6, Ziffern sind, also ganze Zahlen 0,1, ..., 9, sind die 
c; im allgemeinen keine Ziffern mehr. 

Die praktische Durchfiihrung dieser «geordneten» Multiplikation mit Hilfe 
eines Schiebezettels wird als bekannt vorausgesetzt. 

Als abgekiirztes Produkt k-ter Ordnung bezeichnet man die Partialsumme 
bis zum Glied c,, 10* einschlieBlich: : 

(ab) =D cel0t. 
mMrn 
Im Hinblick auf die Verwendung von Rechenmaschinen ist es vorteilhaft, die 
Multiplikation reeller Zahlen unter einem etwas allgemeineren Gesichtspunkt zu 
betrachten. Sei X irgendeine positiv ganzzahlige Potenz von 10, zum Beispiel 
101°, dann kénnen wir unsere reellen Zahlen a, b in der Form 


—co 2 —%7 F 
Get eds. (be BEX! 
Y Ss 


schreiben, wo die A, und B;, jetzt «Ziffern» sind, die zur Zahl X gehoren, also 
ganze Zahlen, die den Ungleichungen 
eee ee 1 OS By Se Xx — 1 


genugen. 
Man kann nun bei passender Wahl von X mit einer Rechenmaschine die Aus- 


driicke 
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bilden und hat auf diese Weise die Méglichkeit, auch mehr abs hundertstellige 
Zahlen ohne groBe Miihe miteinander zu multiplizieren. ; 
Das in diesem allgemeineren Sinn verstandene abgekiirzte Produkt /-ter 
Ordnung der beiden Zahlen a, b bezeichnen wir mit t 
i - 
[ab],= YC, X', ' 
T+S 
wo C, die oben angegebene Bedeutung hat. 

Fiir den Unterschied zwischen dem genauen Produkt a 6 und dem so gewon- 
nenen abgekiirzten Produkt [ab], gelten zwei Abschatzungen, aus denen sich 
fiir den Fall X= 10 die bekannten Fehlerschranken fiir das gewohnliche abge-— 
kiirzte Produkt ergeben. 

Der bei der Bildung des abgekiirzten Produktes [a b], vernachlassigte Rest 
von @ 6 ist gleich 


ab—([a bj, = CL Beta a. Ay 5 137 Se Semel b re B,) Xl-1 
te (Ay By. ns Aggy Biopa eS an Ape 


Der erste Klammerausdruck der rechten Seite enthalt y + s — / 4+- 2 Summanden. 
Wir setzen / = p+ gq (p, g ganz) und behalten auch in den anderen Klammeraus- 
drucken nur die Zahlen A, A"), ..., 4, _,4, und By B, 4, «+3 By o¢4-g oe Alle 
anderen Zahlen A,, B,; ersetzen wir zwecks Abschatzung unseres Unterschiedes 
durch ihren Hoéchstwert X — 1. Ferner setzen wir 


r—st+p S=t+q 
Ae AP id B= eB 
t= j=s 


Dann wird 


ab — [ab], S (A + B) (X — 1) X14 [A+ B+ (X — 1)] (X — 1) X 
+[4+ B+ 2(X — 1)] (X — 1) X88 +>. 


1 1 \ 
=(AF BY(X-YR(14 E+ yet -- +] 
2 2 2 3 fe |S 
+ (Xi xX(14+ 5 4yr4 ) 
2 (X — 1) Xb 1X) + Oe) 
Sak 2) rirararey toe ae =TH x) 
Sy Ws 1B ORM cosas kena ro yee rey 


hae 1)" 
In Worten: 
Bildet man aus den beiden reellen Zahlen 
oo =o: 
a= J A,X* und b= JB, X* 
t=P i=s 
das abgektirzte Produkt /-ter Ordnung [a b],, so unterscheidet sich dieses von 
dem genauen Produkt a b um héchstens (A + B+ 1) X!. Dabei ist A die Summe 
der ersten s— p+ 1A, und B die Summe der ersten r—g+4 1 B,. Die Zahlen 
Aiea: Bo entsprechen einander bei Bildung der /-ten Quersumme. 
Eine besondere Erwaéhnung verdient noch der Fall, daB einer der beiden 
Faktoren, etwa b, nur endlich viele von Null verschiedene Ziffern B, hat: 


t 
b= 5B, xX. 
Ss 


es 
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Dann behalten wir in dem abzuschatzenden Reihenrest alle diese B,;, deren 

Summe wir mit B bezeichnen, bei und ersetzen die A; ausnahmslos durch X — 1. 

“Wir erhalten dann: 

bab —|laby=B (xX — 1) (Xh14 X24 +.) = B(X-— Vor sane Bx 
1 — (1/X) 

In Worten: Die Differenz zwischen dem genauen und dem abgekiirzten Produkt 

7-ter Ordnung der beiden Zahlen 


—oo t 
a= > AgxX* und. b= By X* 
Y Ss 


* * ; ——e nD t 
ist kleiner oder héchstens gleich B X', wo B= J» B, ist. 


Ss 

Der Unterschied zwischen genauem und abgekiirztem Produkt wird natiirlich 
am kleinsten fiir X= 10. A, B und B sind dann Ziffernsummen im gewohnlichen 
Sinne des Wortes. Um nun beide Annehmlichkeiten, namlich die volle Ausntitzung 
der Rechenmaschine und die kleinen Fehlerschranken des gewodhnlichen abge- 
kiirzten Produktes zu genieBen, kann man ein mit der Rechenmaschine, etwa 
unter Zugrundelegung zehnstelliger A,, B,, berechnetes abgekiirztes Produkt 
nachtraglich durch eine kleine erginzende Rechnung mit Schiebezettel in ein ge- 
wohnliches abgekiirztes Produkt verwandeln, dessen Fehler dann wesentlich 
kleiner ausfallen wird. 


(Eingegangen: 9. 3. 1950.) 


EEUU EEE EI ENEEE IEEE 


Buchbesprechungen — Book Reviews — Notices bibliographiques 


i 


Theoretische Mechanik. Von G.Hamev (Springer-Verlag, Berlin 1949). 
796 S., 161 Abb.; DM 63.-. 


Ein ganz ausgezeichnetes Buch, das sich von den iiblichen Lehrbiichern der 
Mechanik durch seinen straffen Aufbau aus einem Minimum an Axiomen unter- 
scheidet. Wie der Verfasser im Vorwort ausfiihrt, war sein erstes Ziel, die ganze 
Mechanik als eine einheitliche Wissenschaft und nicht in der herk6mmlichen 
Weise nach Punktmechanik und Mechanik der Kontinua getrennt erscheinen 
zu lassen. Zweitens sollte die auf LAGRANGE zuriickgehende Methode eines ein- 
heitlichen Aufbaus aus dem Prinzip der virtuellen Arbeiten, dem d’ Alembert- 
schen Prinzip und dem Lagrangeschen Befreiungsprinzip konsequent durchge- 
fiihrt werden. 

Es ist in der Tat erstaunlich, mit welcher Folgerichtigkeit und Eleganz aus 
diesen drei Prinzipien die Mechanik des Massenpunktes, des starren K6rpers und 
des Systems mit endlich vielen Freiheitsgraden entwickelt wird, sodann die- 
jenige der Faden, Seile, Stabe, Haute, Schalen sowie Platten und schlieBlich 
die Hauptsatze der dreidimensionalen Elastizitats-, Hydro- und Aeromechanik. 
Besonders wertvoll ist dabei die ausfiihrliche Diskussion der verwendeten Vor- 
aussetzungen, die dem Leser durch Gegenbeispiele nahegebracht werden. DaB 
im ersten Teil «Aufbau der theoretischen Mechanik», mit seinen 525 Seiten fiir 
die meisten Gebiete nur die Grundlagen dargestellt werden konnten, ist selbst- 
verstandlich. Beim starren Kérper (Kreisel) und beim Punkthaufen (Mehr- 
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kérperproblem) sté&t der Verfasser indessen erheblich weiter vor, und zudem_ 
gibt er eine umfassende Darstellung der unter der Bezeichnung «analytische 
- Mechanik» bekannten Methoden. | 

Der zweite Teil «Aufgaben und Probleme der theoretischen Mechanik» um-_ 
faBt auf 263 Seiten eine reichhaltige Problemsammlung, vom Verfasser im Vor- — 
wort mit der Bemerkung angekiindigt: «Sicher ist Denken wichtiger als Rech- } 
nen. Aber rechnen mu8 man auch kénnen, und man lernt es nicht, wenn man 
daran vorbeigeht. Nichts gefahrlicher als die Meinung, ein Problem sei gelost, 
wenn man die Methode kennt, es zu lésen.» 

Die Lektiire dieses Werkes ist jedermann warmstens zu empfehlen, der mit 
den Elementen vertraut ist. Einmal vermittelt es in selten klarer Weise den 
AnschluB8 an die anspruchsvolleren Werke. Dariiber hinaus ist es aber wie kein 
zweites dazu angetan, Begriffe und Satze zu klaren, die in den elementaren Lehr- 
biichern — teils aus guten Griinden, oft aber ohne Notwendigkeit — mangelhaft 
behandelt werden. 

Dem zuletzt genannten Zwecke dient vor allem das einleitende Kapitel, das 
im wesentlichen eine Axiomatik enthalt. Die Wahl der Axiome ist dabei natur- 
gemaB bis zu einem gewissen Grade willkiirlich; immerhin scheint aber dem 
Referenten das Reaktionsprinzip zu kurz zu kommen, wenn es nicht als Axiom 
formuliert, sondern nur ab und zu beilaufig postuliert und gelegentlich gar als — 
Folge anderer Prinzipien dargestellt wird. Wenn es etwa auf S. 25 und 511 in 
einem Sonderfall auf das Newtonsche Gravitationsgesetz zuriickgefitihrt werden 
kann, das zuvor fiir den ruhenden Zentralkoérper aus den Keplerschen Gesetzen 
erschlossen wird, so spricht dies nur fiir das Gravitationsgesetz in seiner Ver- 
allgemeinerung auf beliebig bewegte Massen, nicht aber gegen das Reaktions- 
prinzip als allgemeines Axiom. Und wenn das Reaktionsprinzip in einem einzigen 
Falle (S. 83) scheinbar nicht gilt, so liegt dies nur an seiner unzuldssigen Erwei- 
terung unter Ausschaltung der zwischen zwei Spiegeln reflektierten Photonen. 

Beilaufig sei noch erwahnt, daB in Aufgabe 76 (S. 609) die Tischkante durch 
eine Offnung in der Platte ersetzt werden sollte, um die betrachtete Bewegung 
mit dem Impulssatz in horizontaler Richtung in Einklang zu bringen. H. Ziegler. 


Darstellende Geometrie. Von F. ReuttTer (Verlag G. Braun, Karlsruhe 
1947 und 1949), Bd. I: 140 S., 144 Abb.; DM 7.50, —- Bd. II: 216 S., 210 Abb.; 
DM 12.50. 


Das vorliegende Werk ist aus einer autographierten Vorlesung hervorgegan - 
gen und will in erster Linie mithelfen, die Schwierigkeiten der Nachkriegszeit 
auf dem deutschen Biichermarkt zu iiberbriicken. Der Aufbau dieses Lehrkurses 
weicht etwas vom Rahmen der bekannten Lehrbiicher iiber darstellende Geo- 
metrie ab; da sich der Verfasser vorwiegend an Ingenieure und Architekten wen- 
det, 1aBt er den Anwendungen in Technik und Architektur etwas breiteren Raum. 

Band I: Orthogonale Zweitafelprojektion. In einem einleitenden Kapitel gibt 
der Verfasser zunachst einen Uberblick iiber die verschiedenen Projektionsarten 
und ihre Invarianten. Nach diesen Vorbereitungen wird dann in der iiblichen 
Weise das Verfahren der orthogonalen Zweitafelprojektion aufgebaut. Besondere 
Kapitel sind der Darstellung von Rotations- und Schraubenflachen sowie den 
Durchdringungen gewidmet. Auf S. 66 findet sich im Text eine kleine Unklarheit; 
im Zusammenhang mit dem ebenen Schnitt eines allgemeinen Kreiskegels ist dort 
von den Dandelinschen Kugeln die Rede, welche ja nur beim geraden Kreiskegel 
vorhanden sind. Die Definition der Ordnung einer Kurve auf S. 92 ist auch fiir 
ein Lehrbuch der darstellenden Geometrie zu unprazis. Darnach heiBt beispiels- 
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-weise eine ebene Kurve von 2. Ordnung, wenn sie mit einer Geraden héchstens 
zwei Schnittpunkte gemeinsam hat. 
; Band IT: Kotierte Normalprojektion, orthogonale Axonometrie und Zentral- 
perspektive. Im ersten Teil gibt der Verfasser einen Uberblick tiber die kotierte 
‘Normalprojektion und ihre Anwendungen in Kartographie und Bauwesen. Bei 
der Behandlung der orthogonalen Axonometrie im zweiten Teil wird der be- 
kannte Weg iiber das Spurendreieck gewahlt. Der dritte Teil ist der Zentralper- 
spektive gewidmet. Die beiden Methoden der freien und der gebundenen Per- 
spektive werden dabei ausfiihrlich beschrieben. Die zahlreichen Figuren erlauben 
dem Leser ein miiheloses Eindringen in die Konstruktionsverfahren der Per- 
spektive. 
Didaktisch wertvoll ist der Anhang des zweiten Bandchens, in welchem der 
Verfasser an Hand einiger Beispiele eine Gegeniiberstellung der verschiedenen 
Abbildungsverfahren der darstellenden Geometrie versucht. M. Jeger 


Electromagnetic Fields, Theory and Application, Vol. 1: Mapping of 
Fields. By ERNST WEBER (John Wiley & Sons, New York 1950). 590 pp.; $10.—. 


In the preface to this book, the first of two volumes on Electromagnetic Fields, 
Professor WEBER explains that the division of material has been made on the 
basis of assigning to Volume I, all static field problems, and to Volume II the 
treatment of dynamic fields. 

Judging from the scope and arrangement of the first volume, this appears to 

_be a very satisfactory method of division for a work designed primarily for grad- 
uate students. The avowed purpose of the book is to give, for field distributions 
derived from potential functions, ‘‘a rather comprehensive survey of the methods 
of analysis and of results obtained with them”’, and it can be said immediately 
that this object has, in very substantial measure, been achieved. Individual 
preference or prejudice will naturally influence opinion as to the appropriate 
variation of emphasis on the sub-divisions of this very wide subject. The author’s 
bias in the matter is possibly illustrated by the fact that some 250 pages out of 
a total of under 600 are devoted to the study of two-dimensional and three- 
dimensional analytic solutions. This represents an exceptionally thorough treat- 
ment of the more difficult potential problems susceptible of solution by mathe- 
matical analysis. It should however be emphasized that the book has clearly not 
been written exclusively for the advanced specialist. Under the title Fields of 
Simple Geometries are developed solutions to the electric and magnetic field 
problems ordinarily dealt with in general text-books. The method of images also 
receives unusually detailed attention. The treatment is clear and concise and 
is characterized by the interjection of explanatory or cautionary notes on many 
points on which misconception or ambiguity commonly exist and by comments 
as to the scope or special value of the original papers or treatises to which reference 
is made. This last is a most valuable feature of the book; references are given 
on a very extensive scale and are judiciously chosen and assessed. 

Experimental and graphical methods of field plotting and numerical methods 
are discussed, but in less detail, and in the reviewer's opinion constitute a rather 
less satisfactory section. It is, of course, true that only a brief outline of these 
individual topics is possible within the compass of a book reviewing the whole 
field but dealing mainly with analysis. Having regard, however, to the important 
part played, for example, by relaxation methods in the practical computation 
of field distributions, a more extensive development of this subject, in basic prin- 
ciple, might have been given. In this section and in that on experimental mapping 
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methods the references to the literature are perhaps a little less complete thar 
in other portions of the book. 
There is an excellent selection of problems associated with each chapter. Fol- 
lowing the main text there are six appendices, two dealing with symbols and 
units, and containing a table of conversion factors from rationalized M. K. Sz 
units in terms of which the book is written, three dealing with basic mathematical 
fundamentals— Vector Analysis, Bessel Functions, and Legendre Functions—and | 
| 


one giving a general bibliography. ; 
The standard of production of the book is in every way excellent. It is a_ 
valuable new work of reference on field problems. James Greig” 


Ausgleichsrechnung. Von V. Happacu (B. G. Teubner Verlagsgesellschaft, 
Leipzig 1950). 104 S., 26 Abb.; $1.50. 


Das kleine Buch gliedert sich in zwei Teile: A. Fehlertheorie, Seite 6—39;— 
B. Fehlerausgleichsrechnung, Seite 40—98. 

Der erste Teil gibt eine einfache Fehlertheorie. Die vorziiglich ausgewahlten 
Beispiele bieten eine ausgezeichnete Wegleitung, und sie werden auch den mit 
der Methode noch wenig vertrauten Wissenschafter, der mit gemessenen Groen 
zu tun hat, dazu veranlassen, die Fehlertheorie anzuwenden. 

Der zweite Teil, der die Hauptformen der Ausgleichsrechnung vorfiihrt, ist 
elementar behandelt. Nach der Ansicht des Referenten wiirde das Buch betracht- 
lich gewinnen, wenn es von der Differenzierung einer Funktion von Funktionen 
Gebrauch machen wiirde. Dadurch wiirde zum Beispiel auf Seite 42 die Herlei- 
tung der Normalgleichungen aus der Bedingung [vv] = Minimum sehr viel ele- 
ganter ausfallen. Ferner ware es zu empfehlen, in aller Kiirze die Bestimmung 
der Gewichte der Unbekannten mit Hilfe der sogenannten Gewichtskoeffizienten 
vorzufiihren. Die im Buche vorgeftihrte Methode, die Gewichtskoeffizienten der 
Unbekannten durch Umstellen der Unbekannten und mehrmalige Neuauflésung 
der Normalgleichungen zu bestimmen, bedeutet eine betrachtliche Mehrarbeit. 
Die notige Theorie kénnte auf wenigen Seiten untergebracht werden. Sie wiirde 
dann auch erlauben, den mittleren Fehler einer Funktion der Unbekannten oder 
der MeBgréBen zu bestimmen. Ein Sachregister wiirde die Verwendung als Nach- 
schlagewerk sehr erleichtern. 

Von diesen Bemerkungen abgesehen, bietet das vorliegende Buch auf knap- 
pem Raum eine korrekte Einfiihrung in die Ausgleichsrechnung und vor allem 
in die Fehlertheorie. Es kann daher vor allem Physikern warm empfohlen werden. 


F. Baeschlin 


Vektoranalysis. Von S. VALENTINER (Walter de Gruyter, Berlin 1950). 
138 S., 19 Abb.; DM 2.40. 


Es ist sehr zu begriiBen, daB die seit einiger Zeit vergriffenen Géschen-Band- 
chen uber Mathematik nun nach und nach wieder neu gedruckt werden. 

Die vorliegende 7. Auflage des Bandchens iiber Vektoranalysis weicht nicht 
wesentlich von den vorhergehenden ab. Es vermittelt nach einer Einfiihrung in 
die Vektoralgebra die Grundlagen der Vektoranalysis (Rechnen mit dem /-Ope- 
rator, Integralsdtze von Gauss, STOKES und GREEN) und bringt dazu einige aus- 
gewahlte Anwendungen auf die Potentialtheorie, die Hydrodynamik und die 
Elektrodynamik. Die zweckmaBig zusammengestellte Formelsammlung der Vek- 
toranalysis im Anhang wird dem Leser sehr willkommen sein. M. Jeger 
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Vorlesungen iiber Differential und Integralrechnung 


Band II: Differentialrechnung auf dem Gebiete mehrerer Variablen 
484 Seiten mit 55 Fig. Ganzleinen Fr. 67.-, brosch, Fr. 63.— (1950) 


Inhalt: Unendliche Mengen — Funktionen auf Mengen — Unendliche Folgen und 
Reihen — Erganzungen zur Differentialrechnung — Anwendungen der Differential- 
rechnung auf die Analysis — Numerische Rechenmethoden — Bogenlange — Ebene 
Kurven — Raumkurven und Flachen. 

In erster Linie enthalt der zweite Band eine vertiefte Behandlung der grund- 
legenden Begriffe der Analysis, ferner eine verhaltnismaBig ausfiihrliche Behand- 
lung der numerischen Rechenmethoden, die sonst etwas stiefmiitterlich erledigt 
werden. Auch in diesem Band findet man eine groBe Anzahl von neuen Aufgaben, 
die sich sowohl auf die praktischen als auch auf die theoretischen Teile des Lehr- 
stoffes beziehen. 

Der dritte und letzte Band befindet sich in Vorbereitung. 
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Die zweidimensionale Laplace-Transformation 
Eine Einfiihrung in ihre Anwendung zur Losung von Randwertproblemen. nebst 
} Tabellen von aie gion . 
a P > ftxen Pat 
Dierricu Weuurar: _ Gustav Dortscu 
Dr. rer. nat. in Heidelberg f Professor in Freiburg i. Br. 


260 Seiten mit 17 Fig. Ganzleinen Fr. 48.-, brosch. Fr. 39.— 


"Inhalt: I. Teil, von G. Doetsch und D. Voelker. ‘Die grundlegenden Eigenschaften 
der &?-Transformation — Partielle Differentialgleichungen erster und zweiter Ord- 


“nung mit zwei unabhangigen Variablen und Systeme von solchen — Partielle Diffe- 
’ rentialgleichungen mit drei unabhangigen Variablen — Funktionalrelationen und © 
- Reihenentwicklungen. 


I. Teil, von D. Voelker. Tabellen von icchpepenieene von allgemeinen Opera- 
tionen und von speziellen Funktionen — Funktionsliste. 


~~ Das Buch, das zum weitaus groBten Teil aus bisher unpubliziertem Material 4 
' besteht, bringt zum erstenmal eine zusammenfassende Darstellung der zvwei- 


dimensionalen Laplace-Transformationen und wendet sie im I.Teil auf eine Fille 


_ von Problemen an, die ftir Physik und Technik wichtig sind, wahrend der II. Teil 


etwa 400 Abbilder von allgemeinen Operationen und 870 Transformationen von 
speziellen Funktionen enthalt. : 


Handbuch der Laplace-Transformation 


von Gustav Doertscnu, Professor in Freiburg i. Br. 


Erster Band: 
Theorie der Laplace-Transformation 


581 Seiten mit 40 Fig. Ganzleinen Fr. 78.-, brosch. Fr. 74.— 
Inhalt: I, Teil: Allgemeine analytische und funktionentheoretische Eigenschaften 


‘der Laplace-Transformation. 


II. Teil: Die komplexe Unokesentele Das Partialintegral der Laplace-Trans- 
formation — Die Parsevalsche Gleichung — Darstellungsproblem — Weitere Umkehr- 
formeln. 

III. Teil: Die, Cesaroschen arithmetischen Mittel der Dgplaceieanaitenntions 


_ IV.Teil: Die Laplace-Transformation der ganzen Funktionen vom Exponential-_ 
typus — Die zweiseitige Laplace-Transformation von analytischen Funktionen — Die — 


Laplace-Transformation der Funktionen der Klasse L?. . 
V. Teil: Abelsche und Taubersche Satze fiir die Laplace-Transformation und fiir 
das komplexe Umkehrintegral. 


Dieses Werk stellt das bisher vorliegende Material ttber die Laplace-Transfor- 


_ mation, vermehrt um zahlreiche bisher nicht publizierte Ergebnisse des Verfassers, 
in einer Weise dar, die sowohl die Bediirfnisse des reinen Mathematikers wie die des 


Ingenieurs beriicksichtigt. Der zweite Band wird die Anwendungen behandeln. . 


VERLAG BIRKHAUSER BASEL 
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